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Zusammenfassung 



In der vorliegenden Arbeit diskutiere ich die hard spectator interaction Amplitude 
von B — > 7T7T zur nachstfiihrenden Ordnung in QCD (d.h. 0(a 2 s )). Dieser spezielle 
Teil der Amplitude, dessen fiihrende Ordnung bei 0(a s ) beginnt, ist im Rahmen der 
QCD Faktorisierung definiert. QCD Faktorisierung ermoglicht, in fuhrender Ord- 
nung in einer Entwicklung in AqcdA 7 ^ die kurz- und die langreichweitige Physik 
zu trennen, wobei die kurzreichweitige Physik in einer storungstheoretischen En- 
twicklung in a s berechnet werden kann. Gegeniiber anderen Teilen der Amplitude 
erfahren hard spectator interactions formal eine Verstarkung durch die zusatzlich 
zur mf,-Skala hinzutretende hartkollineare Skala ^/Aqcd™^, die zu einem grofieren 
numerischen Wert von a s fiihrt. 

Aus rechentechnischer Sicht liegen die hauptsachlichen Herausforderungen dieser 
Arbeit in der Tatsache begriindet, dass die Feynmanintegrale, mit denen wir es zu 
tun haben, bis zu fiinf aufiere Beine haben und drei unabhangige Skalenverhaltnisse 
enthalten. Diese Feynmanintegrale miissen in Potenzen in Aqcd/^6 entwickelt wer- 
den. Ich werde integration by parts Identitaten vorstellen, mit denen die Anzahl der 
Masterintegrale reduziert werden kann. Ebenso werde ich diskutieren, wie man mit 
Differenzialgleichungsmethoden die Entwicklung der Masterintegrale in Aqcd/^6 
erhalt. Im Anhang ist eine konkrete Implement ierung der integration by parts Iden- 
titaten fur ein Computeralgebrasystem vorhanden. 

SchlieBlich diskutiere ich numerische Sachverhalte, wie die Abhangigkeit der Am- 
plituden von der Renormierungsskala und die Grofie der Verzweigungsverhaltnisse. 
Es wird sich herausstellen das die nachstfiihrende Ordnung der hard spectator in- 
teractions wichtig jedoch klein genug ist, so dass die Giiltigkeit der Storungstheorie 
bestehen bleibt. 
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Abstract 



In the present thesis I discuss the hard spectator interaction amplitude in B — > txtx 
at NLO i.e. at 0{a 2 s ). This special part of the amplitude, whose LO starts at 0(a s ), 
is defined in the framework of QCD factorization. QCD factorization allows to sep- 
arate the short- and the long-distance physics in leading power in an expansion in 
Aq CD /m fe , where the short- distance physics can be calculated in a perturbative ex- 
pansion in a s . Compared to other parts of the amplitude hard spectator interactions 
are formally enhanced by the hard collinear scale a/Aqcd^&, which occurs next to 
the m^-scale and leads to an enhancement of a s . 

From a technical point of view the main challenges of this calculation are due 
to the fact that we have to deal with Feynman integrals that come with up to five 
external legs and with three independent ratios of scales. These Feynman integrals 
have to be expanded in powers of Aqcd/ 'n^b- I will discuss integration by parts iden- 
tities to reduce the number of master integrals and differential equations techniques 
to get their power expansions. A concrete implementation of integration by parts 
identities in a computer algebra system is given in the appendix. 

Finally I discuss numerical issues like scale dependence of the amplitudes and 
branching ratios. It will turn out that the NLO contributions of the hard spectator 
interactions are important but small enough for perturbation theory to be valid. 
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Chapter 1 
Introduction 



The present situation of particle physics is the following. On the one hand we have 
got an extremely successful standard model that describes physics up to energy 
scales current accelerators are able to reach. On the other hand it has limitations 
and problems, e.g. the arbitrariness of the standard model parameters, the fact that 
the Higgs particle has not yet been found, the stabilisation of the Higgs mass under 
loop corrections (fine tuning problem) or the question why the electroweak scale is so 
much lower than the Planck scale (hierarchy problem). So most particle physicists 
expect new physics to show up at energy scales that are beyond the range of present 
accelerators but will be reached by future colliders. Within the next year LHC at 
CERN will start running and in the following years will collect data from proton 
proton collisions at a centre of mass energy of 14 TeV. This will allow us to obtain 
information about new physics by producing not yet observed particles directly. On 
the other hand physics beyond the standard model can be discovered by precision 
measurements of low energy quantities which are influenced by new physics particles 
because of quantum effects. The currently running experiments BaBar and Belle 
and after the start of LHC also LHCb are dedicated to examine decays of 5-mesons, 
where new physics is expected to be seen in CP asymmetries. 

However in order to find new physics by indirect search, some parameters of the 
standard model have to be determined more precisely. To this end LHCb will make 
an important contribution. Above all the Wolfenstein parameters p and fj [U Ej that 
occur in the parametrisation of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
and determine its complex phase, which leads to CP asymmetry in the standard 
model, are up to now only very imprecisely known [3]: 

p = 0.182±85Sf f} = 0.3321H. (1.1) 

These parameters, which influence weak interactions of quarks, can be determined 
with higher accuracy by 5-meson decays. In order to reduce their large uncertainties 
on the experimental side better statistics is needed, which is expected to be improved 
in the next few years, and on the theoretical side we have to get hadronic physics 
under control. This is due to the fact that weak interactions of quarks, from which 
p and fj are measured, are always spoiled by non-perturbative strong interactions, 
because quarks are bound in hadronic states like mesons. Hadronic physics, however, 
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is governed by the energy scale Aq CD , where QCD cannot be handled perturbatively. 

There are several advantages in observing 5-meson decays. One of them is due 
to the production of 5-mesons itself: There exists an extremely clean source to 
produce B- mesons: The resonance T(4S), a bound state of a bb pair, has a mass 
that is only slightly larger than twice the mass of the 5-meson and decays nearly 
completely into BB pairs. This resonance is used at the 5-factories BaBar and 
Belle. 

Another advantage is the possibility to obtain clean information about the com- 
plex phase of the CKM matrix by measurement of quantum mechanical oscillations 
in the B — B system. The lifetime of 5-mesons, which is about 1.5 ps, is large 
enough to observe those oscillations in the detector [H [3] . By measuring the time 
dependent CP violation it is possible to obtain the CKM angles a, (3 and 7, which 
determine p and fj, with small hadronic uncertainties (see e.g. chapter 1 of [6]). The 
"golden channel" B — > J/tpKs, where the dependence on hadronic quantities is 
strongly suppressed by small CKM parameters, allows a quite precise determination 
of sin(2/5) = 0.687 ±0.032 [7]. In the same way the decay B — > nn could be used for 
a precise determination of the angle a. However other than in the "golden channel" 
in the case of B — > nn hadronic physics plays a subdominant but non-negligible 
role. 

At this point another convenient property of the 5-meson comes into play. The 
mass of the 6-quark introduces a hard scale, at which a s is small enough to make 
perturbation theory possible. However the bound state of the fe-quark in the B- 
meson is dominated by physics of the soft scale Aqcd ; where perturbation theory 
breaks down. While inclusive decays can be handled in the framework of operator 
product expansion, for exclusive decays, which the present thesis deals with, the 
framework of QCD factorization has been proposed [HIE]- This framework makes it 
possible to disentangle the soft and hard physics at leading power in an expansion in 
AqcD/ m fe- Decay amplitudes are then obtained in perturbative expansions, which 
come with hadronic parameters that have to be determined in experiment or by non 
perturbative methods like QCD sum rules or lattice QCD. 

Whereas the a s corrections for the transition matrix elements of B — > titi have 
been calculated in [10J, the present thesis deals with the 0(a 2 s ) contribution of a 
specific part of the amplitude. This part, which consists of the hard spectator 
interaction Feynman diagrams, will be defined in the next section. There I will also 
argue, that it is reasonable to consider the hard spectator interactions separately. 
The calculation of the rest of the 0(a 2 s ) corrections has been partly performed by 
Guido Bell in his PhD thesis [TT| IT2"] . There the complete imaginary part and a 
preliminary result of the real part of the amplitude is given. My calculation of the 
hard spectator scattering amplitude is not the first one as it has been calculated 
recently by [TBI III] . It is however the first pure QCD calculation, whereas [131 EH] 
used the framework of soft-collinear effective theory (SCET) [To], [TBI E7J an effective 
theory, where the expansion in Aqcd/^6 is performed at the level of the Lagrangian 
rather than of Feynman integrals. It is the main result of this thesis to confirm the 
results of [131 E] an d to show by explicit calculation that pure QCD and SCET 
lead to the same result in this special case. 
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From a technical point of view the calculation in this thesis consists of the eval- 
uation of about 60 one-loop Feynman diagrams. The challenges of this task are due 
to the fact that these diagrams come with up to five external legs and three inde- 
pendent ratios of scales. In order to reduce the number of master integrals and to 
perform power expansions of the Feynman integrals, integration by parts methods 
and differential equation techniques will prove appropriate tools. Most parts of the 
calculation will be performed by a computer algebra system, whereas the algorithms 
and the necessary steps to obtain input results for the programs will be discussed in 
detail. As I did not obtain the completed 0(a 2 s ) corrections, the phenomenological 
part of this thesis is restricted to the reproduction of the branching ratios numeri- 
cally obtained in [T3]. Other observables like CP asymmetries are not improved by 
my partial result alone. 

This thesis is organised as follows: In chapter [2] I start with an introduction 
to QCD factorization and define in this framework the hard spectator scattering 
amplitude. After defining my notations I demonstrate the calculation of the LO 
of the hard spectator interactions and end the chapter by explaining the technical 
details of the integration by parts methods and differential equation techniques. 

Chapter [3] is the most technical of all. There all of the Feynman diagrams that 
contribute are listed and their evaluation is discussed in detail. Furthermore NLO 
corrections to the wave functions and evanescent operators occurring at this order 
are dealt with. 

After presenting the complete analytical results and the numerical analysis in 
chapter H] I end up with the conclusions. 
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Chapter 2 



Preliminaries 



2.1 Hard spectator interactions and QCD factor- 
ization 

Though the decay of the 5-meson is caused by weak interactions, strong interactions 
play a dominant role. It is however not possible to handle the QCD effects completely 
perturbatively. This is due to the energy scales that are contained in the 5-meson: 
Whereas a s at the mass of the 6-quark is a small parameter, the bound state of 
the quarks leads to an energy scale of (9(Aq CD ) which spoils perturbation theory. 
The idea of QCD factorization [HI E] is to separate these scales. At leading power in 
Aqcd/^6 we obtain the amplitude for B — » tttt in the following form: 



Two different types of quantities enter this formula. On the one hand the hadronic 
physics is contained in the form factor F B ~* n and the wave functions 4>bi and 
which will be defined in the next section more precisely. These quantities contain the 
information about the bound states of the mesons. They have to be determined by 
non-perturbative methods like QCD sum rules or lattice calculations. Alternatively, 
because they are at least partly process independent, they might be extracted in 
the future from experiment. On the other hand the hard scattering kernels T 1 
and T 11 contain the physics of the hard scale 0(rrib) and the hard collinear scale 
0( a/^Aq CD ) and can be calculated perturbatively. 
Here I would like to make two remarks to (12. II) : 

First I want to note that (12.11) is only valid in leading power in the expansion 
in Aqcd/ 7 ^?!- Higher orders in this expansion lead to endpoint singularities i.e. 
the integrals over the variables x, y and £ diverge at the endpoints. This leads to a 
mixture of the physics of the soft scale into T 1 and T 11 and spoils QCD factorization. 
However there are corrections that are formally of subleading power but numerically 




(2.1) 



6 



2. Preliminaries 





Figure 2.2: Hard spectator interactions at 0(a s ). This is the LO of T 



enhanced and cannot be handled within the framework of QCD factorization. They 
have to be estimated in the numerical analysis. 

The second remark concerns the separation of the hard scattering kernel into 
T 1 and T n . The Feynman diagrams that contribute to B — ► can be distributed 
into two different classes. The class of diagrams where there is no hard interaction 
of the spectator quark (fig. 12.11) contributes to T 1 . The hard spectator scattering 
diagrams, which are shown in LO in a s in fig. 12. 2\ contribute to T u . Through 
the soft momentum I of the constituent quark of the B-meson the hard collinear 
scale ^/AqcD m 6 comes into play. This leads to the fact that in contrast to T l , 
which is completely governed by the scale m&, T n has to be evaluated at the hard- 
collinear scale. This leads to an enhancement of a s and makes the NLO corrections 
(i.e. 0(a 2 s )) of the hard spectator interaction diagrams more important. These a 2 s 
corrections of T 11 are the topic of the present thesis. Hard spectator scattering 
corrections to the penguin diagram (third diagram of fig. 12. ip are beyond the scope 
of this thesis. The cancellation of the dependence on the renormalisation scale 
of this class of diagrams is completely independent of the "tree amplitude" i.e. 
the diagrams of fig. 12.21 and higher order a s corrections. For phenomenological 
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applications, however, they should be taken into account. 

2.2 Notation and basic formulas 
2.2.1 Kinematics 

For the process B — > tctt we will assign the momenta p and q to the pions which 
fulfil the condition 

p 2 ,q 2 = 0. (2.2) 

This is the leading power approximation in Aqcd/ 777 ^ as we count the mass of the 
pion as O(Aqcd)- Let us define two Lorentz vectors n + , n_ by: 

< = (1,0, 0,1), nt = (1,0, 0,-1). (2.3) 

In the rest frame of the decaying meson p can be defined to be in the direction of 
n + and q to be in the direction of n_. Light cone coordinates for the Lorentz vector 
z^ are defined by: 



z° + z 3 z° - z 3 



So one can decompose z^ 1 into: 



z x = (0,z\z 2 } 0) (2.4) 



= ^P q » + llljf + z » (2.5) 
p • q p ■ q 

such that 

^± • V = z ± ' 1 — 0- (2.6) 
We denote the mass of the S-meson with and the mass of the 6-quark with 
mi,. The difference tub — rrib = O(Aqcd) such that we cannot distinguish those 
masses in leading power. However setting 

m b = m B (2.7) 

in Feynman integrals might lead to additional infrared divergences. So we have to 
perform the integral before we can make the substitution (12.71) unless we are sure 
that we do not produce infrared divergences. If we calculate Feynman integrals, it 
is convenient to set 

m B = 1 (2.8) 

such that p • q = \. The dependence on m# can be reconstructed by giving the 
correct mass dimension to the physical quantities. 

2.2.2 Colour factors 

In our calculations we will use the following three colour factors, which arise from 
the SU(3) algebra: 

Cn = \, °p = ^^- and °g = N c , (2.9) 
where N r = 3 is the number of colours. 
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2.2.3 Meson wave functions 

The pion light cone distribution amplitude 0^ is defined by 

(*toMz) a [..tf(0)f,\0)* =o = ^(M^^dxe^M*)- (2-10) 
The ellipsis [. . .] stands for the Wilson line 

[z,0]=Pexp(J dtig s z ■ A(zt)\ , (2.11) 

which makes (12.10j) gauge invariant. For the definition of the S-meson wave function 
(j>Bi we need the special kinematics of the process. Following [9] let us define 

*g(z, PB ) = (0Mz)[. ■ .]6 a (0)|B(p fl )> = / j^e-^fdpjs). (2.12) 



(2tt) 4 

In the calculation of matrix elements we get terms like: 

' tr(A(/)<M0) = / T^Z / d*ze a -*tr(A(l)* B {z)). (2.13) 



(2tt) 4 v v " J (2tt) 4 j 

We will only consider the case that the dependence of the amplitude A on I is like 
this: 

A(l) = A(2l-p) (2.14) 

In this case we can use the 5-meson wave function on the light cone which is given 
by 0: 



(2.15) 

2 ,zx=0 



(0\q a (z)[. . .}b p (0)\B(p B )) 

Jo 



where ^ 

/ d£$ B1 (0 = l and f d^ B2 (O = 0- (2.16) 
Jo Jo 

It is now straight forward to write down the momentum projector of the 5-meson: 
^tr(A(2/.p)¥(l)) 



tr(^ + m B ) 75 / ^($ m (0 + ^ + $B 2 (0)^(^l) (2-17) 
•/ o 



1 

At this point we give the following definitions 

m B 



/f&nCO (2-18) 
A n = / *ln»^si(0- (2-19) 



2.2 Notation and basic formulas 
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2.2.4 Effective weak Hamiltonian 

The effective weak Hamiltonian which leads to B — > nn is given by [TH] : 



G 



p=u,c 



i=3...6 



+ h.C, 



(2.20) 



where = V* d V pb and 



01 


= (dp) v _ A (pb)v_ A , 


(2.21) 




= (diPj)v-A(pjbi)v-A, 


(2.22) 


o 3 


= (db) V -A^2(qq)v-A, 
q 


(2.23) 


o A 


= {dib j )v-A S ^ J {.qjqi)v-A, 

q 


(2.24) 


o 5 


= (db) v _ A ^2(qq) v+Aj 

q 


(2.25) 


o 6 


i 


(2.26) 


o 8g 




(2.27) 



Explicit expressions for the short-distance coefficients C; can be obtained from [18 
The decay amplitude of B — >• irn is given by 



A{B -> TXTl) = (7TTT\H ca \B). 

For later convenience we define 



A(B 



A(B 



7T7T , 



A{B 



7T7T 



(2.28) 



(2.29) 



where A 1 {A 11 ) belongs to the first (second) term of (12.1 1) . Because A 1 and A 11 
contain different hadronic quantities, the renormalisation scale dependence of both 
of them has to vanish separately. So we can set their scales to different values /i 1 
and jtz 11 . As in A 1 there occurs only the mass scale m& we can set /i 1 = m&. 111 A 11 
there occurs also the hard-collinear scale ^/AqcD m b- As we will see this scale is an 
appropriate choice for /i 11 . 

In order to separate the QCD effects from the weak physics we write the matrix 
elements of the effective weak Hamiltonian in the following factorised form [TO] : 



\B) 



(2.30) 



p=u,c 
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Figure 2.3: Annihilation topology. The gluon vertex that is marked by the cross 
can alternatively be attached to other crosses. 

where 

T p = ai5 pu {ub) v _ A <g> (du) V -A 
+a 2 5 pu (db) v _ A (g) (uu)v-a 

+a 3 ^2(db) V - A ® (qq)v-A 
q 

+«4 ^2(qb)v-A ® (dq)v-A 
q 

+a 5 ^2(db) V - A <8> (gg)y+A 

+<^(-2)(#) 5 „p ® (Jg) 5+ P- (2.31) 
q 

Note that in contrast to [TU] the electroweak corrections to the effective weak Hamil- 
tonian are not included in the above equations as in the case of B — > tttt they can 
be safely neglected. T p ann stands for the contributions of the annihilation topologies, 
which are shown in fig. 12.31 These contributions do not occur in leading power 
and cannot be calculated in a model independent way in the framework of QCD 
factorization. For the exact definition of 71 ann I refer to |10j . The matrix ele- 
ments of the operators j\ ® j 2 are defined to be {ixit\ji <8> J2I-B) = ( 7I "|ji|-B)( 7r |j2|0) 
or (7r| j'2 |-B)(7r|ji|0) corresponding to the flavour structure of the 7r-mesons. The 
penguin contractions that are shown in the third diagram of fig. 12.11 and the contri- 
butions of the operators O^-Ogg are by definition contained in the amplitudes 0,3-0%. 
As we take in the present thesis only the "tree amplitude" (fig. I2.2p into account, 
we only calculate a 2 s corrections to the amplitudes a\ and a 2 . 

The decay amplitudes of B — > tttt can be written in terms of a« as follows [TO] : 

-A(B°^<K + ir-) = [A^ + ^K + r^H, 
-V2A(B~ -> n-n°) = X' u ( ai + a 2 )A nw 

A(B°^tt°tt ) = [-\' u a 2 + \' p (a p 4 + ry 6 )] A nn (2.32) 

where 

A^ = i%{m 2 B -ml)f^U 



2.3 Hard spectator interactions at LO 
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and 



2mi 



x m h {ji){mjji) + m d (ji))' 
For the LO and NLO results of the a, I refer to [TOl. 



(2.33) 



The annihilation contributions are parametrised in the following form [TO] : 



A (B° - 
A (B° 



7T 7T 



7T-7T ) 

► ttV) 



where 



[A^ + ^ + A'J (63 + 264)] B„ 




(2.34) 



(2.35) 



The parameters 6j can be further parametrised by the Wilson coefficients occurring 
in (12.201) and purely hadronic quantities: 



b, 
64 



c 3 Ai + c 5 (4 + 4) + 7V C C 6 4 



(2.36) 



where the quantities A 1 ^ are given by [TO] : 



4 



7TCt, 



18 X A - 4 + — 



2rfXi 



4 


127ta s rl (2X 2 A - X A ) . 



(2.37) 



Here r£ is defined as in (I2.33f) and Xa parametrises an integral that is divergent 
because of endpoint singularities. In section 14.4.21 1 will give an estimate of Xa for 
numerical calculations. 



2.3 Hard spectator interactions at LO 

The leading order of the hard spectator interactions which start at 0(a s ) is shown 
in fig. 12. 2[ The hard spectator scattering kernel T n , which does not depend on 
the wave functions, can be obtained by calculating the transition matrix element 
between free external quarks, to which we assign the momenta shown in fig. 12.21 
The variables x, x = 1 — x, y, y = 1 — y are the arguments of T , which arise from 
the projection on the pion wave function (I2.10p . In the sense of power counting we 
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count all components of I of (9(Aqcd) ; while the components of p and q are (D(ra&) 
or exactly zero. We define the following quantities 

£ = Lit 

~ p-q 

6 = (2.38) 
p-q 

We will see that in the end the dependence on 9 vanishes in leading power such that 
we can use ( I2.17p . 

We consider the three cases B° — > ir + iT~, B° — > 7r°7r° and B~ — > 7r~7r°. In the 
case, that the external quarks come with the flavour content of I? — > tt + tt~ , the LO 
hard spectator amplitude for the effective operator O2 reads: 

(i) ( m - 



^spcct. (B 



7T 7T 



(d(xp)u(xp) u(yq)d(yq)\02\d{l)b(p + q — I)) 



spect. 



1 



Ana s C F N c — —d(l)Yd(xp) u(xp)Y(l ~ l$)b{p + q — I) 



d(yq) ^ 



x£m 

( 



y yy 



z7m7v (l-7sMi/g), (2-39) 



where the quark antiquark states in the input and output channels of the matrix 
element form colour singlets. The subscript "spect." means that only diagrams with 
a hard spectator interaction are taken into account. The amplitude of 0\ vanishes 
to this order in a s . In the case of B° — > 7r°7r° we get the tree amplitude from 
the matrix element of 0\. The case B~ — ► tt^tt does not need to be considered 
separately, because from isospin symmetry follows [T9], [10] : 

y/2A(B~ -> 7T"7T ) = A(B° -> 7r + 7r) + A(B° -> 7T 7T ). (2.40) 

On the other hand the full amplitude is the convolution of T 11 with the wave 
functions, given by (12.11) . To extract T n from (12.391) we need the wave functions 
with the same external states we have used in (12.39}) . i.e. we have to calculate the 
matrix elements ( I2.10p and ( 12.121) . where the pion or B- meson states are replaced 
by free external quark states. To the order O{a s ) we get 

<t> ( Xp(y') = y"^(^g)e-^'(«(z/g)rf(i/g)l4^)<(0)|0) z - i2±=0 

= 27rN c 6(y' - y)d p (yq)u a (yq) 
^lapi x ') = 2nN c 5(x' - x)u p {xp)d ct {xp) (2.41) 

Ap{l'-) = Jdz + e a '- z+ {0\de(z)%(0y\d(l)b(p + q-l)) z -, zx = 
= 2irN c 5(l'~ - r)d p (l)b a (p + q - I) 

By using 

4pL = / ^"^Lw^^^ftfO^fei/^)^ ( 2 - 42 ) 



2.3 Hard spectator interactions at LO 
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we finally obtain: 



C 



T II(1) (x,y,/ ) a < a/W7 = 4:ira s \^ / [7^(1 - 7 5 )] a , 



2 ' 7 L » V / 0/ Ja'7 



(2yr) 3 iV c 2 

(2.43) 



t7m7^ (1 - 75j 



y yy 



It should be noted that only the first summand of the above equation contributes 
after performing the Dirac trace in four dimensions. The second summand is evanes- 
cent. This will be important, when we will calculate the NLO corrections of the wave 
functions (see section [3731) . 

If we plug the hadronic wave functions defined by (12.101) and (I2.17P into (I2.42p i.e. 
we calculate the matrix element (12.391) between meson states instead of free quark 
states, we get for the LO amplitude □: 

4pL = - %f ^lf/ ^ f dxd V^ $Bl(fl6r(s)6r(y)-^- (2-44) 

Following ( 12. ip and the conventions of p3] we write our amplitude in the form: 

A spcct .i = -im 2 B dxdyd£Tj\x,y,£)f B $Bi(Ofn(pn(x)fn(pn(y)- (2.45) 
Jo 

where in the case of B — > 7r + 7r~ we define 

^spect.l = (G^spect. 
Aspect. 2 = 

spcct. 

(2.46) 

and in the case B — > 7r°7r° we define 

Aspect. 1 (C^l)spect. 

Aspect. 2 — (C^spcct.- (2.47) 

Because we use the NDR-scheme which preserves Fierz transformations for 0\ and 
C>2, Tf has the same form for both decay channels. From (12.441) and (I2.45P we get: 

T ii(i) A C F 1_ 

1 ~ s 4N^xym 2 B 
T 2 II(1) = . (2.48) 

According to pU] the contribution of (12.441) to a± and a 2 (see (12.311) ) is given by: 

C 2 C F 7ia s 

0*1.11 — ^ tljnr 

N c 

dCpiras 

02,u = ^2 H ™ [2A§) 



1 At this point I want to apologise to the reader because of a mismatch in my notation: A spoc t. 
is used for the matrix elements of the effective Operators Oi between both free external quarks 
and hadronic meson states. It should become clear from the context what is actually meant. 
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where 

= J4f f f *bi(0 /' ^r(*) (2.50) 
m BJ+ Jo ? >/ o Jo y 

and the label II in ( 12. 49ft denotes the contribution to .A 11 as defined in ( 12.29ft . 



2.4 Calculation techniques for Feynman integrals 
2.4.1 Integration by parts method 

Integration by parts (IBP) identities were introduced in j20j [21] . An algorithm to 
reduce Feynman integrals by IBP-identities to master integrals is very well described 
in [22]- So I will only show the basic principles. Because the topic of my thesis is a 
one-loop calculation I will restrict to the one-loop case, the generalisation to multi 
loop is straight forward. 

The most general form0 of a one-loop Feynman integral is 

d d k (k ■ p h ) ni . . . (k ■ p k ) ni 



(2n) d [(k + p h ) 2 - Mir ■■■[(* + P.) 2 - M?p 

1 f d d k s^.-.sf 1 



k- Pli +M( 



mi 

2 



lll-u 



(27r) d D™ 1 . . . D^Df 1 . . . Df 



(2.51) 



where n x ■ ■ ■ n u m 1 . . . m t , m x . . . m u > 0, j\ . . . , j u h . . . , i t , h, . . . , i u G {1, ...,n} 
and pi, . . . ,p n are the momenta which appear in the internal propagator lines. With- 
out loss of generality we can assume that there is no k 2 in the numerator as we can 
make the replacement 

k 2 = Di + Ml - p? - 2k ■ p h . (2.52) 

Because of our special kinematics we have only three linearly independent momenta 
p, q, I, so all of the momenta p±, . . . ,p n are linear combinations of p, q, I. This will 
simplify the reduction of the Feynman integrals. We will define 

B = { Pl ,...,p k } (2.53) 

to be a basis of span{pi, . . . ,p n } where k < 3 and pi, . . . ,pk G {pi, . . . ,p n }- 

Following [22] we can reduce (I2.5ip by performing algebraic transformations on 
the integrands, which are defined in the following three rules: 

Rule 1. Consider the case that there exist {ci, . . . , q} such that 

t t 
^c j p i] = and J^c,- = 1. (2.54) 
j=i 3=1 



2 Tensor integrals which contain expressions like fc M , k^k" , . . . in the numerator can be reduced 
to scalar integrals as described in (23] 
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Now we can make the following simplification (I 6 {1 . . . , t}): 
k-p it 



D™ 1 ■■■ D? 



lA 

2 

s Eft 



Pi 



i - Ci 



+ 



D mt 



(2.55) 



and the scalar product k ■ p^ has disappeared from the numerator. If ( 12. 54ft cannot 
be fulfilled we use the identity 

'Tf-M 2 1+ pl-pl) + k- Pil 



k-p k = -(A - D x 



(2.56) 



to reduce our set of integrals further. This identity does not reduce the total number 
of scalar products in the numerator but the number of different scalar products. 

Rule 2. For scalar products of the form k ■ pj. we make the replacement 



k ■ PL 



D j 



Ml 
bf 



(2.57) 



(2.581 



Rule 3. Now consider the case that our integrand is of the form 

k-Pk 

D™ 1 . . . D^Df 1 . . . Df* 

where pk ^ {Ph, ■ ■ ■ ,Pi t jPi 1 T ■ ■ iP\ }■ I n that case we use the following rule: Choose 
a set bi C {p^ l ,...,pi } which is a basis of span{pj , . . . ,P~ iu }- Choose b 2 C 
{p iiy . . . ,p it } such that b = b x U b 2 forms a basis of spanjjo^, . . . ,Pi t ,pi , ■ ■ ■ ,Pj u }- 
Complete b to a basis of span{p 1; . . . ,p n } by adding elements of {pi, . . . ,p n } to 
b. Then write pk as a linear combination of this new basis and apply (if possible) 
(12351) . (12~56|) or (12~57D respectively. 

For the following identities which are called integration by parts or IBP identi- 
ties we will use the fact that in dimensional regularisation an integral over a total 
derivative with respect to the loop momentum vanishes. Using the definitions of 
( 12.511) we get two further rules: 

Rule 4. 

d d k d k^sl 1 ...s? 







{2n) d dk^ D[ ni . . . DT bf 1 . . . Df 



(d 



2r) 



d d k 



. . . s 



in 



(2n) d Df 1 ...D^Df 1 ...Df 



a=l 



mi 



Dt 



d d k 

{2Tt) d 
u 

-E 

a=l 



{Ml 



. . s 



in. 



k- Pia s? 



,ni 



PI 



11 . . . D ma+1 . . . DT 



D™ 1 . . . D\ 



■ ■ DT 



d d k 



k-P~ in 



{27r) d bf 1 . . . bf^ 1 ...bi 



(2.59) 
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where s = X)i=i n « an( ^ r — Y^i=i m «- 

Another identity is: 
Rule 5. 

ci d A; 9 p^ 1 . . . s|" 







(27r) d dfc/ 1 D[ ni . . . D™* Df 1 . . . Df» 

^ ™6Pa ' P * y ( 27 r)d £)™i . . . £>m tpm x ^ 

d d k si 1 ... s™< (A; • p a + p ib ■ Pa) 



6=1 

t 



6=1 "* 



(2Ti) d D™ 1 ... D™ b+1 . . . DT'Df 1 ...Df 
d d k sT ... s? 1 



6=1 



m bPa ■ P ib J ^ ^ ^ ^ ^ ^ ^ ^ D ^ b+1 ^ ^ ^ D ^ (2.60) 



where p a e B. 



For the IBP identities (12. 59ft and (12.601) we have used the translation invariance 
of the dimensional regularised integral. We get another class of identities if we use 
the invariance under Lorentz transformations. From equation (2.9) of [23] we get 

/ d d d d 



dpi "^dpl ™"dp? t ^dpl 

d d d d 

^ ilV dp- dp- P~^ v Qp? PiuH QpL> 

d d d d \ 

Pnu dpj ~ Vnii dp^ x + • • • + Pnu dpj t ~ Pn "M ) x 

d d k s? 1 . . . s™ ; 



(2ir) d D?\..D™Df 1 ...Df™' 
We choose Pi,Pj G B. By multiplying of (12.611) with pfpj we get 
Rule 6. 



(2.61) 



r I 



b 1 . . . b a . . . b l 



D™ 1 . . . Dl nt D™ 1 . . . D™ 



a=l 1 • ' ' a • ■ ■ ^\ • • • ii, 

S 1 1 ... S,' 



™a(k ■ PiP- a -Pj-k- p jPla ■ Pi ) -— 



o=l 



D mi . . . D mt D mi . . . D™^ 1 . . . jyhu 



(2.62) 



An implementation in Mathematica of the IBP identities can be found in ap- 
pendix [A] 
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2.4.2 Calculation of Feynman diagrams with differential equa- 
tions 

In this section I will discuss the extraction of subleading powers of Feynman integrals 
with the method of differential equations [25l l26l [24] . This method will prove to be 
easy to implement in a computer algebra system. The idea to obtain the analytic 
expansion of Feynman integrals by tracing them back to differential equations has 
first been proposed in (25]. This method, which is demonstrated in (25] by the one- 
loop two-point integral and in |26j by the two-loop sunrise diagram, uses differential 
equations with respect to the small or large parameter, in which the integral has to 
be expanded. 

In contrast to [251 EE] I will discuss the case that setting the small parameter 
to zero gives rise to new divergences. In this case the initial condition is not given 
by the differential equation itself and also cannot be obtained by calculation of the 
simpler integral that is defined by setting the expansion parameter to zero. It is 
not possible to give a general proof, but it seems to be a rule, that one needs the 
leading power as a "boundary condition" , which can be calculated by the method of 
regions [23 [28], [221 [30] . The subleading powers can be obtained from the differential 
equation. In the present section I will discuss which conditions the differential 
equation has to fulfil in order for this to work. 

Description of the method 

We start with a (scalar) integral of the form 



where the propagators are of the form Di = (k + pf) 2 — mf. We assume that there 
is only one mass hierarchy, i.e. there are two masses m <C M such that all of the 
momenta and masses pi and m, are of 0(m) or of O(M). We expand ( 12. 63ft in jj 
by replacing all small momenta and masses by Pi — > Xpi and expand in A. After the 
expansion the bookkeeping parameter A can be set to 1. 

We obtain a differential equation for I by differentiating the integrand in ( 12.631 ) 
with respect to A. This gives rise to new Feynman integrals with propagators of 
the form -=j and scalar products k ■ pi in the numerator. Those Feynman integrals, 

i 

however, can be reduced to the original integral and to simpler integrals (i.e. integrals 
that contain less propagators in the denominator) by using integration by parts 
identities. 

Finally we obtain for (12.63P a differential equation of the form 



where h(X) contains only rational functions of A and g(X) can be expressed by 
Feynman integrals with a reduced number of propagators. It is easy to see that h 
and g are unique if and only if I and the integrals contained in g are master integrals 




(2.63) 



d 

dX 



I(X) — h(X)I(X) + g(X) 



(2.64) 
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with respect to IBP-identities, i.e. they cannot be reduced to simpler integrals by 
IBP-identities. If J(A) is divergent in e = /, h and g have to be expanded in e: 

i = £V 

i 

/i = J^e* 

i 

9 = J>e f . ( 2 - 65 ) 

i 

Plugging (12.651) into (I2.64p gives a system of differential equations for Jj, similar to 
( 12.641) . In the next paragraph we will consider an example for this case. 

First let us assume that h(X) and g(X) have the following asymptotic behaviour 
in A: 

h(X) = M 0) + Xh® + . . . 

g(X) = ^A^WQnA) (2.66) 

3 

i.e. h starts at A , and we allow that g starts at a negative power of A. We count 
In A as C(A°) so the g^> may depend on In A. This dependence, however, has to be 
such that 

limA# (j) (lnA) = 0. (2.67) 
The condition (I2.67P is fulfilled, if the are of the form of a finite sum 



^a n ln"A. (2.68) 



n=no 



The limit m — > oo however can spoil the expansion (12361) . E.g. e" lnA = I so the 
condition ( 12.671) is not fulfilled, which is due to the fact that we must not change 
the order of the limits A — * and m — > oo. 

Further we assume that also /(A) starts at A 

/(A) = J (0) (In A) + A/« (In A) + . . . (2.69) 

and plug this into (12.641) such that we obtain an equation which gives 1^ recursively: 

A*J« = J* dX'X'^ 1 h^I^' 1 ^ (In X') + g^XlnX'^J . (2.70) 

I want to stress that, because h starts at O(X ), (12.701) is a recurrence relation, i.e. 
i"w does not mix into /W if j > z. As the integral is only well defined if i > 1, we 
need the leading power as "boundary condition" and ( 12.701) will give us all the 
higher powers in A. It is easy to implement ( 12.701) in a computer algebra system, 
because we just need the integration of polynomials and finite powers of logarithms. 
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A modification is needed if h starts at A 1 i.e. 



h = -^h^ + .... (2.71) 
A 



By replacing / = X n I we obtain the differential equation 



i 1= f^ + h) I + X n g (2.72) 



d\ \\ 
which is similar to f)2.64p and leads to 

/A /i+n-l \ 

dX'X'^ 71 " 1 ^ h®I<^ 1 -n (In A') + g {i - 1] (In A') J , (2.73) 

which is valid for i > 1 — n. So, if I starts at 0(\~ n ), the subleading powers result 
from the leading power. 

Examples 

We start with a pedagogic example: 
Example 2.4.1. 

f d d k 1 

where Ac 1. The exact expression for this integral is given by: 

lnA(l + A + A 2 + ...). (2.75) 



(An) 2 1 - A (4tt) 2 

We see that / diverges for A — > 0. As described e.g. in [30] we can obtain the leading 
power by expanding the integrand in the regions k ~ vA and k ~ 1. This leads in 
the first region to 

r(l + 6)fi + l-lnA^ (2.76) 



(27r) d k 2 (k 2 - A) (4vr) 2 ^ 
and in the second region to 

d d k 1 z 



(2-K) d k\k 2 - 1) (47r) 2 - e v V e 
such that we finally obtain 



r(l + e) - + 1 (2.77) 



/(°)(lnA) = ^lnA. (2.78) 

This is the result we obtain from the leading power of (12.751) . We write the derivative 
of / with respect to A in the following form: 
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We obtained the right hand side of (12.791) by decomposing 4rl into partial fractions. 
Of course this decomposition is not unique which is due to the fact that / itself 
is not a master integral but can be further simplified by partial fractioning. From 
(EZD and (12341) we get: 

h = — *— = 1 + A + A 2 + ... 
1 — A 

9 = (4^A|T^A) = (4^( A " + 1 + A + -) < 2 ' 80 > 

such that the coefficients in the expansion in A according to ( 12.661) do not depend 
on the power label (k): 

h {k) = 1 and g {k) = (2.81) 

(4"7T ) l 



We obtain for the recurrence relation (12.701) : 



Using the initial value ( 12.781) it is easy to prove by induction 

I ik) (In A) = In A \/k > 0. (2.83) 

[An) 2 

This result coincides with (I2.75p . 

The first nontrivial example, we want to consider, is the following three-point 
integral: 

Example 2.4.2. 

/= / d d k 1 

J (2n) d k 2 (k + un^ + l) 2 (k + n + + n^) 2 ' 1 ' ; 

Here n + and n_ are collinear Lorentz vectors, which fulfil = n 2 _ = and 
n+ -n_ = |, u is a real number between and 1 and / is a Lorentz vector with I 2 = 
and P <C 1. Furthermore we define 

£ = 2Z-n+ and = 21 ■ n_. (2.85) 

We expand / in /, so we make the replacement / — > XI and differentiate / with respect 
to A. The integral is not divergent in e such that we obtain a differential equation 
of the form f!2.64j) where the Taylor series of h(\) starts at A as in (I2.66p . In g(X) 
only two-point integrals occur, which are easy to calculate. I do not want to give 
the explicit expressions for h and g because they are complicated, their exact form 
is not needed to understand this example and they can be handled by a computer 
algebra system. Because the leading power of / is of C(A°), ( I2.70P gives all of the 
subleading powers. 
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We obtain the leading power as follows: First we have to identify the regions, 
which contribute at leading power. If we decompose k into 

F = 2k ■ n+nt + 2k ■ n_n\ + k^ (2.86) 

we note that the only regions, which remain at leading power, are the hard region 
k^ L ~ 1 and the hard-collinear region 

k ■ n + ~ 1 
k ■ n_ ~ A 

~ V\. (2.87) 

The soft region k^ ~ A leads at leading power to a scaleless integral, which vanishes 
in dimensional regularisation. In the hard region we expand the integrand to 

(2.88) 



k 2 (k + unJ) 2 {k + n + + n_) 2 



By introducing a convenient Feynman parametrisation we obtain for the (4 — 2e)- 
dimensional integral over (12. 88ft : 

T(l + e) exp(zTre)- ( Ml^l _ I ln 2 (1 _ u) \ . ( 2 .89) 



(4vr) 2 - v ' v-^v— / u 
In the hard-collinear region we expand the integrand to 

1 



k 2 (k + un_ + 6n + ) 2 {2k ■ n+ + 1) ' 

The integral over (12.901) gives: 

7 , \ / . 1 /- ln(l — u) . . 1, 9 . , 

T(l + e) exp(zTre)- ^ '- + 2Li 2 (u) + - ln 2 (l - u) 



(2.90) 



(4tt) 2 - v ' ~>—^"-> u \ e ' 1 2 

+ lnuln(l -w) + ln(l -u)ln0 1. (2.91) 



So adding f)2.89p and f)2.9ip together we get the leading power of (12. 84ft : 

j(°) = — i_ -(2Li 2 (u)+lnuln(l-«)+ln(l -u)ln0). (2.92) 
(47r) 2 m 

By plugging (12~92|) into (l2~70l) we obtain J at 0(A): 
/« = 

i 1L/ „ , ln(l-u)ln0 ln(l-u)lnu , , 2Li 2 (wl 
01 -2 + lnu + — + — + ln£ + 



(4tt) 2 u 



u u u 



Inu | ln(l — w) | ln(l — tt)ln0 | ln(l — u)lnu 



1 — u u u u 

lnf 2Li 2 (w) 



1 — It M 

Now we want to consider the following four-point integral 



(2.93) 
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Example 2.4.3. 

f d d k 1 

~ J (2vr) d k 2 {k + n„) 2 {k + I - n + ) 2 {k + I - un + ) 2 ' ^ > 

where we used the same variables, which were introduced in (12.841) . This example 
is very special, because in this case our method will allow us to obtain not only the 
sub leading but also the leading power in I. J is divergent in e such that we obtain 
after the expansion (I2.65P a system of differential equations of the following form: 

d r 

-77-1-1 = M-l + 9-1 

-^Io = h I + h x I-i + go- (2.95) 

It turns out that in our example h takes the simple form 

2 + 2e 



(2.96) 



A 

such that analogously to (I2.72p we can transform (I2.95P into 

>/-,) = AV, 

^(A%) = -2A/_ 1 + A 2 9o. (2.97) 

This system of differential equations can easily be integrated to: 

1 r x 

(») _ 1 / (*-!) 



J« = J^^ A dA'r +1 (-2J« + ^) (2.98) 

where the superscript (i) denotes the order in A as in (I2.66P and (12.691) . Both 7_x 
and Jo start at (9(A _1 ). Because (12.981) is valid for i > —1, it gives us the leading 
power expression, which reads: 

(47r) 2 " e u£ \e 1 - u J 

where £ = 21 ■ n + as in the example above. The exact expression for (12.941) can be 
obtained from [21]. Thereby (12.991) can be tested. 



A simplification for the calculation of the leading power 

In the last paragraph I want to return to Example 12.4.21 I will show how we can 
use differential equations to prove that the integral ( 12. 84ft depends in leading power 
only on the soft kinematical variable 9 = 21 ■ n_ and not on £ = 21 ■ n + . We 
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need derivatives of the integral with respect to £ and 9, which we have to express 
through derivatives with respect to Z M . These derivatives can be applied directly 
to the integrand, whose dependence on is obvious. We start from the following 
identities: 

a d _ d_ T ,,_d_ T 

n+ dl» ~ 89 + hp 

"*-4rJ = il + O-IL 1 ( 2 - 10 °) 



which lead to 



6 k = 5(-K+f"t + n^ 

where we have set I 2 = in (12.1011) . Using (12.1011) we can show that in leading; 
power (12.841) depends only on 9 and not on £. So we can simplify the calculation of 
the leading power by making the replacement Z M — > 9n+. The proof goes as follows: 
From ( 12.691) we see that the statement "J*- -* does not depend on £" is equivalent to 

&I(£\,6\) = 0{\). (2.102) 

Using the first equation of (12.1011) we get 

^/(£A, 9X) = 0(A)I(£A, 9X) + 0(A). (2.103) 

Because we know (e.g. from power counting) that /(£A, 9X) starts at A , (I2.102p is 
proven. 
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Chapter 3 

Calculation of the NLO 



3.1 Notation 
3.1.1 Dirac structure 

In this thesis I used the NDR scheme [32] such that 75-matrices are anticommuting. 
We get for the matrix elements of 0\ and 2 Dirac structures of the following type: 

(0i,2) = gi(/)rigi(xp) q 2 (xp)T 2 b(p + q-l) q~z{yq)Tzqi{m) (3.1) 

where qi are u- or <i-quarks. To avoid to specify the flavour, which depends on the 
decay mode, I introduce for (13. ip the following short notation: 



The equations of motion lead to: 



i®r 2 c 


sr 3 








®r 2 


®r 3 


= 


Fit 


®r 2 


®r 3 


= 


rid 


vpT 2 


®r 3 


= 


-/+ 


m b ) 


®r 3 


= 


rid 


>r 2 <s 




= 




>r 2 <s 


5T 3 ^ 


= 



(3.2) 



(3.3) 



3.1.2 Imaginary part of the propagators 



Unless otherwise stated propagators in Feynman integrals always contain an term 
+17] where 77 > and we take the limit 77 — * after the integration. For example an 
integral of the form 

d d k 1 



(2vr) d k 2 (k + p) 



2 



is just an abbreviation for 

d d k 



(2vr) d (k 2 + ir))({k + p) 2 + irj)' 
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If the propagator does not contain the integration momentum quadratically, the irj 
will always be given explicitly. 



3.2 Evaluation of the Feynman diagrams 

The diagrams that contribute to T± at NLO are listed in fig. 13. 11 - 13. 71 Fig. l3.1l shows 
the gluon self energy. This is a subclass of diagrams which factorizes separately. 
After adding the counter term for the gluon propagator the contribution to T* 1 ^ is: 



7-n(2) 



AN^xy 



C 



N 



16 



mi 



lnf 



16 



hire 



80 



.'5u 2 5, , 5, 31 



(3.4) 



For (13.41) we have set the number of active quark flavours to rif = 5 and set the 
mass of the u-, d-, s- and c-quark to zero. 

For the rest of this section I will consider only those diagrams, for which the 
calculation of Feynman integrals is not straight forward. I will only show how those 
Feynman integrals can be calculated in leading power. For higher powers I refer to 
the methods shown in section 12.4.21 It is important to note, that though for the 
evaluation of the Feynman integrals arguments depending on power counting have 
been used, all of the integrals occuring in this section have been tested by methods 
that do not depend on power counting. 

The first class of diagrams with non trivial Feynman integrals are the diagrams 
in fig. 13.31 The two diagrams from fig. 13.3( a) read 

alii = -ig 4 s N c C F (C F - ~C G ) 

f d d k 7M (^-/) 7r 7 ,(l- 75 )0 7 ^(^ + ^+y^-/)7^(l- 75 )(^-^ 
J (2ix) d k 2 (k-l) 2 {k + xp-l) 2 (k-yq) 2 (k + xp + yq-l) 2 

(3.5) 



aII2 = -ig 4 s N c C 2 F 




d d k 7 „(ff - /) 7 r®7.(l - 7s) ® l^yj- #)7"(1 ~ 7s)(ff + V<1+ - /)7 r 
(2vr) d k 2 (k - l) 2 {k + xp - l) 2 {k - yq) 2 {k + xp + yq - I) 2 



(3.6) 

The denominators of (13.51) and (13.61) are identical up to the substitution y —> y. So 
the Feynman integrals we have to calculate are the same. We can reduce our five- 
point integrals to four-point integrals by expanding the denominator of the integrand 
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Figure 3.3: Diagrams all 
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Figure 3.5: Diagrams alV 
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Figure 3.7: Nonabelian diagrams 
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into partial fractions 



k 2 {k — l) 2 (k + xp — l) 2 {k — yq) 2 {k + xp + yq — 1)* 



y(x - [k 2 (k - l) 2 (k + xp- l) 2 (k - yq) 2 

V 1 

yk 2 {k — l) 2 {k + xp — l) 2 {k + xp + yq — I) 2 
1 

k 2 (k — l) 2 (k — yq) 2 (k + xp + yq — I) 2 

y i 

y (k — l) 2 {k + xp — l) 2 {k — yq) 2 {k + xp + yq — I) 2 



(3.7) 



Only the first two summands of the right hand side of this equation give leading 
power contributions to alii and aII2: 



The third one 



ti = 
tu = 

tin 



1 



k 2 (k - l) 2 (k + xp - l) 2 (k - yq) 2 
1 

k 2 {k — l) 2 {k + xp — l) 2 {k + xp + yq — I) 2 



1 



k 2 {k — \) 2 {k — yq) 2 {k + xp + yq — I) 2 
gives only a leading power contribution in the hard-collinear region 

k ■ p ~ 1 
k ■ q ~ A 



(3.8) 
(3.9) 

(3.10) 



and the soft region 



P ~ A, 



(3.11) 



(3.12) 



where we introduced the counting ~ A and set ms — 1- In both regions the 
leading power of the numerators of alii and aII2 vanishes because of equations of 
motion. The fourth summand of ( 13.71) 



tiv 



1 



(k — l) 2 (k + xp — l) 2 (k — yq) 2 (k + xp + yq — I)' 



(3.13) 



does not give a leading power contribution at all. 

The topologies £j and tu are exactly the denominators of the Feynman integrals 
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of the last four diagrams of fig. 13.31 

1 f d d k 



k 2 {k — l) 2 (k + xp — l) 2 (k + xp + yq — I) 
1 „ . 1 f d d k 



2 



aII6 = ig A s N c C 2 F __ \ x (3.17) 
xy — xt, — yti 



aII4 = - ^c)-— ^ 7 " (3.15) 

^ 7,(l-75)®7"(l-75)(^^+y^- /h^ + ^+^i/- /)7 r 
k 2 {k — l) 2 {k + xp — Z) 2 (/c + xp + yq — I) 2 

aII5 = irfJV^CCfc - ±— ^ / -0^^ - /) 7r (3.16) 

^ 7.(1-75)®7 M (Z/^- ^)l T {xp / +y4- /) 7 "(l- 75 ) 
/c 2 (£; - /) 2 (fc + - /) 2 (/c - yg) 2 

'F~ 

xy 

d d k j E (l- /) 7T< i, 7y (l- 75 )g) 7 ^(l- 75 )(x^ + ^- /)7 r (y^- ^7 M 
(27r) d P(fc-Z) 2 (A; + xp-Z) 2 (£;-yg) 2 

where the Feynman diagrams are given in the same order as they occur in fig. 13.31 
So we have reduced all the Feynman integrals of fig. 13.31 to the topologies ti and tu- 
Regarding tj we need the scalar integral 

Doi{l) s / wr h - (3 - 18) 

By following the procedure of (23] the tensor integrals J m^kHj, J -^ak^k^tj and 



(2w) dn -<' J (27r)° 

J -^^k^k u k T ti can be reduced to (13.181) and to the two-point master integrals that 
are listed in appendix IB. 1[ The leading power of (13.181) gets only contributions from 
the region where k is soft i.e. all components of k are of (9(Aqcd)- In this region we 
can expand the integrand of (13. 18ft to 

(3.19) 



{k 2 + irj)((k — I) 2 + irf){2k ■ p — £ + irj){— 2k ■ q + irf)xy 

Now we can obtain the leading power of (13.181) by integrating (13.191) over all momenta 
k because there is no other region, which gives a leading power contribution, besides 
where k is soft. The integration of (13.191) can be easily performed by using the 
following Feynman parametrisation |33j : 

, TTn+T " (3-20) 
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3. Calculation of the NLO 



Finally we obtain for the leading power of (13.181) 



D 



01 



i r(l +e)(47r/^) e / 2 21n£ + 21n# + 2z7r 



(An) 2 xytf 



- 7T 2 + ln 2 £ + ln^ # + 21n£ln# + 2ni(\n^ + ln( 



(3-21] 



Regarding the topology t n we need the tensor integrals J j^LkHii, J -^ 3: k^k u t II 



(27r) d 



and J -^dk^k u k T tu. This topology can be reduced to two-point master integrals 
and to the four-point master integral 



f d d k 



(3.22) 



In order to calculate this integral we decompose tn into 



1 



y [k 2 (k - l) 2 (k + xp- l) 2 (2k -q + x-8) 

1 

~k 2 (k - l) 2 (k + xp + yq- l) 2 (2k -q + x 



(3.23) 



where only the integral over the first summand of (13. 23f) gives a leading power contri- 
bution. This integration is straight forward if one uses the Feynman parametrisation 

(E2DD. 

Finally we get the leading power of (13.221) : 



D 



on 



r(i 



i)(A-KjJL 



2\e 



(An)' 



x 2 yi 



2 tt 2 

-(lnx + ln£) h ln 2 x + ln 2 £ + 21nxln£ 

6 3 



(3.24) 



Actually it turns out that the sum of the diagrams in fig. 13.31 vanishes in leading 
power. 

The diagrams of fig. 13.41 are straight forward to calculate. It is easy to see that in 
leading power I does not occur within a loop integral. So there are only two linearly 
independent momenta in the Feynman integrals, which, using similar relations like 
(13. 7ft and IBP identities, can be reduced to the master integrals that are listed in 
appendix IB. 1L 

The diagrams of fig. I3.5( a).(b) are easy to calculate, because they contain only 
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three-point integrals. The next two (fig. 13.5( c)) are given by 



alVml = ig A s N c C F (C F - -C G ) 



d d k 
(2tt) c 



7^(^ + x^)7 r ® 



(3.25) 



7 "(l- 75 )(^+^+ j- / +mft ) 7r (g )7y (l- 75 )(^- f + xp' + y^ % 

k 2 (k + xp) 2 {k + xp — l) 2 (k + xp + yq — l) 2 {k 2 + 2k ■ (p + q — I)) 



aIVm2 = -ig 4 s N c C F (C F - -C G ) 



d d k 

W) 



x 



(3.26) 



7"(1 -75)(# + ?f+ ^~ /+m e ,)7 T <8)7 M (^- l+xj + yjfryQ -75) 
A; 2 (/c + xp) 2 {k + xp — l) 2 (k + xp + yq — l) 2 {k 2 + 2k ■ (p + q — I)) 

By expanding the denominator of (13.251) or (I3.26P into partial fractions we get: 

1 

k 2 {k + xp) 2 ik + xp — l) 2 {k + xp + yq — l) 2 {k 2 + 2k ■ {p + q — /)) 



x — x£ — 9 
1 



fc 2 (/i; + xp) 2 (k + xp — l) 2 (k + xp + yq — I) 2 
1 



+ 



yk 2 (k + xp) 2 (k + xp- l) 2 (k 2 + 2k- (p + q-l)) 

V 1 

yk 2 {k + xp) 2 {k + xp + yq — l) 2 (k 2 + 2k ■ {p + q — I)) 
x 1 

x k 2 {k + xp — l) 2 (k + xp + yq — l) 2 {k 2 + 2k ■ (p + q — I)) 
x 1 

x (k + xp) 2 {k + xp — l) 2 (k + xp + yq — l) 2 {k 2 + 2k ■ (p + q — /)) 
where we get leading power contributions only from: 

1 



(3.27) 



tn 
tin 



k 2 {k + xp) 2 (k + xp — l) 2 (k + xp + yq — I) 2 

1 

k 2 {k + xp) 2 {k + xp- l) 2 {k 2 + 2k ■ (p + q - I)) 

1 

(k + xp) 2 {k + xp — l) 2 {k + xp + yq — l) 2 {k 2 + 2k ■ {p + q — I)) 



(3.28) 
(3.29) 
.(3.30) 



The leading power of ( 13.28P can be taken from (13.211) . We obtain the leading power 
of (I3.29P by making the replacement I — » £q. Alternatively we can use ( ]B.18j) in ap- 
pendix IB. 21 and take the leading power afterwards. For (13.301) we obtain the leading 
power by making the replacement I — > 9p. In contrast to (I3.29P this replacement is 
not so obvious and to calculate this integral exactly is very involved. However we 
can start with the value, we obtained by this prescription, and show afterwards by 
solving a partial differential equation that it is correct. We define 



/d d k 
—Tjt in (x,y,p,q, XI) 



(3.31) 
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3. Calculation of the NLO 



and derive two differential equations by deriving I with respect to x and to y. As a 
boundary condition for our differential equations we calculate / at the point x = 
and y = 1. This can be done by decomposing tm into partial fractions and using 
IBP identities. We use the fact that the limits x — > and y — > 1 do not lead to 
extra divergences in e and in A. So we can solve our differential equations order by 
order in e and A. Defining 

I = 5^lj*VA* (3.32) 

j,k 

and using the fact that we only need the leading power in A we obtain differential 
equations of the form 

Ir(-i) _ h r(-i) + /, r(-i) + „ 



<9x 

l/j-" = (3-33) 

The coefficients h x0 , h xl , g X j \ h yo and g V j are straight forward to calculate 
using IBP identities and the master integrals are given in appendix IB. 11 It turns 
out that the leading power integral we obtained by the prescription I — > Op fulfils our 
boundary condition for x = and y = 1 as well as the set of differential equations 

The sum of the following two diagrams (fig. 13.5( d)) is 




-tg s l\ c U F {U F — —L>g 



2 w 7 {2Ti) d k 2 (k-l) 2 (k-xp) 2 (k + xp~l) 



2 



^ (k — xp — yq) 2 ^ u ^ u (k — xp — yq) 2 ^^ ) 

As in the previous case the denominator can be decomposed into partial fractions 
and the remaining four-point master integrals can be calculated in leading power by 
the replacement I — > £q. Alternatively we can use the explicit formulas for four-point 
integrals given in [31] and derive the leading power (and higher powers if necessary) 
afterwards. This is what I have done in order to get an independent test of my 
master integrals. 

The last two diagrams of this class (fig. 13.5( e)) are given by 

i r i r 
i- + — 



igtN c C 2 F 



d d k 7^7 T g)7 r ^(l -7 5 ) 
{2n) d k 2 {k + xp) 2 (k + p) 2 (k + I) 2 



^(k + l + yq) 2 ^ ~ ^(k + l + yq) 2 ^ ' (1 " ^ ^ 
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As in the case before we obtain the leading power of the four-point master integrals 
by the replacement I — >• £q or by using the formulas in [31]. The denominator of 
( 13.351) however cannot by decomposed into partial fractions. So we additionally need 
the five-point master integral given in appendix IB. 31 

Regarding the diagrams in fig. 13.61 there do not occur any subtleties we have not 
yet considered in the paragraph above. So we directly switch to the non-abelian 
diagrams in fig. 13.71 The first four (fig. 13. 7( a). (b)) do not lead to any problems 
because they contain only three-point integrals. The diagrams in fig. 13.7( c) are 
given by 



(3.36) 



-tgtN c C F C 



d d k 
x£ J (2vrp 
/-0) 7 "(l- 75 



(g" x (l - k) T + 2g Xr k» - g Ti "{k + 21- 2xp) x ) x + 



k 2 (k + l-p) 2 (k + l-xp) 2 



<8> 



7s)- 



(k + yq) 2 (k + yq) 2 

The four-point integral we have to solve is nearly the same as (12.941) of example 
12.4.31 We can reduce this integral to a solution of a differential equation and get in 
this way every power in AQ CD /m b . 

The last diagrams we will consider are those from fig. 13.7( d). In leading power 
they read: 



(3.37) 




ig*N c C F C G - 

X^ J (27T' d 

7^7^(1 -7s)(# + ^ 



(g^(k-xpY 
^+1)7a _ 



2g Xr k il + g Tii {k + 2xp) x ) x 



k 2 (k + xp-l) 2 ((k + p + q) 2 - 
-{„(% + x$ + yd- /) 7r 7t 



xp' + yd- /) 7l / 



(k + xp + yq 
The scalar master integral 
d d k 



[k + xp + yq — l)' 



;i-7 5 ) 



(2n) d k 2 ik + xp- l) 2 (k + xp + yq - l) 2 {{k + p + q) 2 - 1] 



(3.38) 



can be calculated in leading power by setting 6 = i.e. we make the replacement 
\p — > ^qp. This can be seen as follows: Counting soft momenta as O(X) and hard 
momenta as 0(1) (remember = 1) the regions of space where (I3.38P gives a 
leading power contribution are 

F ~ 1 
fc" ~ A 

K ■ p ~ A fcV ~ k ■ q ~ 1. 
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In these regions occurs only in the combination I - p. So we can make the replace- 
ment / M — > Those people who do not believe these arguments are invited to use 
the exact expression ( IB. 171) for the four-point integral with one massive propagator 
line, which is given in appendix IB. 21 After taking the leading power it can easily be 
seen that we get the same result as by just making the replacement / M — > 

The diagrams which contribute to TJ, 1 are those of fig. 13.31 and fig. 13.41 The 
other diagrams drop out because their colour trace is zero. As in the case of T* 1 
the diagrams of fig. 13.31 cancel each other in leading power. The remaining diagrams 
are easy to calculate because their Feynman integrals are the same as in the case of 



3.3 Wave function contributions 
3.3.1 General remarks 

It has already been demonstrated in section 12.31 how in principle we can extract the 
scattering kernel T 11 of (12. ID from the amplitude if we know the wave functions. T 11 
does not depend on the hadronic physics and on the form of the wave function <p n 
and <pB in particular, so we can get T 11 by calculating the matrix elements of the 
effective operators between free quark states carrying the momenta shown in fig. 12.21 
on page EJ Because we calculate T 11 in NLO we need unlike as in section 12.31 the 
wave functions up to NLO. Let us write the second term of ( 12.1 ft in the following 
formal way: 

Apcct. = 0,r ® fa ® <Pb ® T u . (3.39) 

All of the objects arising in (13.391) have their perturbative series in a s , so (13.391) 
becomes 

<L, = 0f ® # ® <4 0) ® (3.40) 

4pL = € ] ® # ® <4 0) ® + # ® ® <4 0) ® t 11 ^ + 



where the superscript (i) denotes the order in aJE In order to get T ll ^> we have to 
calculate .Aspect.) 4$ an d 4$ f° r our ^ na ^ states. Then T ll< ^> is given by 

# ® # ® ^ 0) ® T n ( 2 ) = (3.41) 
Agl, - # ® # ® <4 0) ® - # ® 0? ) ® 0g) ® - 
0f ® 0f ® 0« ® T n « 

At this point a subtlety occurs. Let us have a closer look to the factorization formula 
(12. ip . By calculating the first order in a s of the partonic form factor F B ^ 7T '( 1 \ which 



1 Please note that the hard spectator scattering kernel starts at 0(a s ). So we call T 11 ^ 1 ' the 
LO and T 11 ^ the NLO. 
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Figure 3.8: Example for diagrams which obviously belong to the form factor 



is defined by free quark states instead of hadronic external states, we see that it can 
be written in the form 



pB-nr,(l) = (O) 3 (O) T 



<1) 

' formfact. 



(3.42) 



But T f ^ mfact is no part of T u . So we have to modify (I3.4ip insofar as we have to 



subtract the right hand side of (13.421) from the right hand side of (13.411) : 



(3.43) 



A {2) 



if 



T ii(i) 



T ii(i) 



,(o) 



r (i) 

" formfact. 



T I(1) 



In (13.431) we did not include the term _F B ^ 7r >( 2 )(g)T I ( ' ) , because it is obviously identical 
with the diagrams where the gluons do not interact with the emitted pion (e.g. those 
of fig. 13.81) . Those diagrams where not considered in the last section. So we do not 
have to consider them here. 

The wave functions for free external quark states are given at LO by (I2.4ip . At 
NLO there exist three possible contractions: The two external quark states can be 
connected by a gluon propagator or one of the external quarks can be connected to 
the eikonal Wilson line of the wave function (fig. 13.91) . The diagrams of fig. 13.9( a). (b) 
and (c) give 0(ot s ) of the "pion wave function for free quarks", i.e. we have replaced 
the pion final state (k(p)\ in (12.101) by the free quark state (q'(xp)q(xp)\. The Fourier 
transformed wave function 0^ {%') is defined analogously to (12.411) . For the diagrams 
in fig. 13.91 (a), (b) and (c) respectively we get: 



J -KOLf3 



X 



]n 2 ia s C F N c 
ln 2 ia s C F N c 
]n 2 ia s C F N c 



d d k 8{x' - x - K) - 5{x' - x) _ 

Qf3[xp) 



(2vr 



k 2 k~ 



d d k S(x - x' - |+) - 5{x - x') 



(2vr) c 



k 2 k~ 



d d k S(x' -x + £) 



(2k) 



k 2 



q(xp)j' r 



xi 



1 



x\ 

1 



7 5 (xp) 



q' a (xp) 



t — X'j 

——j T q'(xp) 
+ xp 



(3.44) 
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b(p + q-l) q(l) 
(d) 



(e) 



(f) 



Figure 3.9: NLO contributions to the meson wave functions. The dashed line stands 
for the eikonal Wilson line which makes the wave functions gauge invariant. 

3.3.2 Evanescent operators 

At NLO the convolution of the wave functions with the tree level kernel T n 'W gives 
rise to new Dirac structures, which, however, can in four dimensions be reduced to 
the tree level Dirac structures. So we obtain the tree level Dirac structures plus 
further evanescent structures, which vanish for d = 4 but give finite contributions if 
they are multiplied by UV-poles. We define our renormalisation scheme such that 
we subtract the UV-poles and these finite parts of the evanescent structures. 

The tree level kernel (I2.43[) contains two Dirac structures where the second one 
is evanescent (after the the projection on the wave functions). We write T 11 ^ in the 
following form (that must not be mixed up with the notation (13.21) ): 

T u ^(x,y,r) = J-y®7"(l - 75) ® (^T ~ (1 - 75) (3.45) 

xi v y yy J 

where the symbol <§) stands for the "wrong contraction" of the Dirac indices i.e. the 
Dirac indices are given by 

[r^r 2 ®r 3 ] a , a/3W7 = r^, 7 r^ (346) 

as in (12.431) . The "right contraction" is defined by the symbol <8> i.e. writing the 
Dirac indices explicitly 

In d = 4 the wrong and the right contraction are related by Fierz transformations. 
It is convenient and commonly used to define the renormalised wave functions in 
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terms of the right contraction, i.e. to define </>^ cn " by renormalising the operator 
g(z)7 /1 75g'(0) instead of q(z)pq'(Q) a . This is why we define our renormalisation 
scheme such that only the UV-finite part of the right contraction operators remains: 
Using the notation of (13.451) - ( 13. 47ft we define the following operators: 

1 2i> 
O (x,y,r) = --—7^(1 -7s)® 7/,(l -75)®— (l-7s) (3.48) 
Ziii cc y 

Ox{x,y,l-) = 7^7^(1-75)®— (1- 75) (3.49) 
rx y 

2 (x,y,r) = J r7 ^(l-7 5 )®Z^M^(i_ 75 ). (3.5O) 
l~x yy 

The matrix elements of these operators are defined analogously to (12.421) : 

(Oi) = J dxdy'dl'- (f) vaa '{x) ^wiv') <fiB~ r Y>(l'~)Oi a > a { 3 >py~ / (x' ,y' (3.51) 

Note that (0\ + O2) is just the convolution of the tree level kernel (I3.45j) with the 
wave functions. Furthermore by using Fierz identities it is easy to prove that we 
have in four dimensions 



(Oo) = (Ci) (3.52) 

(0 2 ) = 0. (3.53) 

So we define the following evanescent operators: 

E 1 = 2 (3.54) 

E 2 = Oi-Oo (3.55) 

If (2- d) 2 \ 
Ez ^ x~y~F (W^W^ 1 - 7s) + 2 7 M (1 - 75) ® 7m(! - 75) J 

(8)2^(1 - 75) (3.56) 

E A = -^ 7r 7 M 7V®7r7A7"(l - 7s) ® ^7,(1 - 7s) (3.57) 



where we have defined £3 and E± for later convenience. Using those operator def- 
initions we define our renormalisation scheme such that we subtract the UV-pole 
of (Oo) and the finite parts of (Ei) i.e. terms of the form -^-(E), where (E) is 
an arbitrary evanescent structure. It is important to note that we do not subtract 
IR-poles, because they depend not only on the operator but also on the external 
states the operator is sandwiched in between. They have to vanish in (13.431) such 
that the hard scattering kernel is finite. Finally we obtain the same result as if we 
had regularised the IR-divergences by small quark and gluon masses because the 
evanescent structures vanish in d — 4. 

In the next step we will calculate the convolution integral of T 11 '^ with the 
NLO wave functions given by (13.441) . i.e. we have to calculate the renormalised 
matrix elements of 0\ + 02 at NLO. 
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3.3.3 Wave function of the emitted pion 

First we consider the renormalisation of the emitted pion wave function: Because the 
contribution of the wave functions <fi^ and (p£ does not change the Dirac structure 
of the operators, we do not need to consider evanescent operators when we calculate 
the diagrams of fig. 13.9( a). (b). So for the emitted pion wave function these diagrams 
give after renormalisation: 



/nren.\(l),(a),(b) 
2 / emitted 



2a. 

47T 



em Mir / V 

- Q- + 21n// IR ) (2 + ln2/ + lnj7)(0 2 )C°) 



(3.58) 

where the LO matrix elements (Oi)^ can be obtained from ( 12.391) . Note that we 
kept the IR-pole times the evanescent matrix element (O2) • This is needed for 
consistency because we also kept similar terms in the QCD-calculation of v4 spect .. 
Furthermore it allows us to show that all IR-divergences vanish. 

The diagram in fig. 13.9( c) mixes different Dirac structures. So we have to include 
evanescent operators in the renormalisation. In the case of the emitted pion wave 
function the operator 0\ does not mix under renormalisation with the evanescent 
operator E\ 03.54p . We obtain for the renormalised matrix element: 



(0[ en -) 



(l),(c) 
emitted 



2a s y]ny 



1 



— + 2 In 

e IR /^IR 



A.ix y 

The matrix element of E\ however has an overlap with 0\ 



(Oi) {0) . 



(3.59) 



IF \ (1) ' (c) 

X-^l /emitted 



^-^ + 21n^X 
euv em A*ir 



-2y\xiy-2y\xiy) (Ek)<® 
y\ny s 



— 4e In y 



<(9i> (0) 



(3.60) 



The renormalisation prescription tells us to subtract the UV-pole and the UV-finite 
part of (Ei). So we obtain after renormalisation: 



/ ciren.\(l),(c) 
/emitted 



Air' 



c f 



_± + 21n^X ) {-2y\ny 
em /im 

+ 4(\ny + ^)(O i r 



2y\ny) {E,) 



(o) 



(3.61) 



Note that the evanescent operator E\ leads to a finite term 4(lny + ^-^■)(0\) 
which we would have missed if we had just dropped the evanescent operators. 



>(o) 
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3.3.4 Wave function of the recoiled pion 

In the next step we consider the NLO contribution of the recoiled pion wave function. 
As in the case of the emitted pion the diagrams fig. 13. 9( a). (b) do not lead to a mixing 
between the operators. Therefore we get: 



, //ycn.\(l),(a),(b) 
^2 /recoiled 



a s 2 1n x + 4x 

~—C F 

Att x 



-— + 21n^ 
em Mir 



(Oi) (0) 



(—4 > 2 > (0) 



(3.62) 



Other than in the case of the emitted pion the operators 0\ and O2 mix the 
spinors of the recoiled pion and the 5-meson. Therefore we have to work in the 
operator basis of Oq and the evanescent operators and define our renormalisation 
scheme such that the finite parts of the matrix elements of the evanescent operators 
vanish. The diagram fig. 13.9( c) contributes to the matrix element of the renormalised 
operator C^ 611 ': 



/£jren.\(l),(c) 
\ / recoiled 



2^f-L-21n^)^(O > (0) . 



47T 



em 



Mm 



x 



(3.63) 



In the case of the evanescent operators we keep the IR-pole: 



/ Diren.\(l),(c) 
\ 1 / recoiled 



/ iriren.\(l),(c) 
\- C/ 2 /recoiled 



la S/1 ( \ n . \ x\yyx , ^ . 

-— C F — + 21n/x m £4 

2 4vr V. e iR / x 

la s „ / 1 , \ irlnx . ^ . 
4 47T V e iR / x 



(0) 



(0) 



(3.64) 
(3.65) 



At the end of the day we obtain a contribution from diagram fig. 



//?-)ren. , rnren. , rnren.\ (l),( c ) 
\ U + ^1 + ^2 /recoiled 



\ a s x In x 

— 1> f 

2 4tt x 



— -21n^X 
e m Mir 



1 

"xl~ 



2jig 



nviv V(o) 
yy 



8<Oi> 



(o) 



(3.66) 



The very complicated but also very explicit form, in which the above equation was 
given, is rather convenient, because on the QCD side the diagrams in fig. 13.5( e) on 
page [281 come with the same Dirac structure and cancel the IR-pole of (13. 66ft . 
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3.3.5 Wave function of the B- meson 



The a s corrections of the wave function of the S-meson are given by the second row 
of fig. 13.91 For the diagrams (d), (e) and (f) respectively they read: 



87r 2 ia s N c C F 

<t>$$\n = 

8n 2 ia s N c C F 



k 2 k~ 



d d k 6(1'- -I- - k~) - 6(1'- - r) 
(2vr) rf 

d d k 6(1'- 



I. 



(3.67) 
b a (p + q-l) 



l~ — k~) — 6(l'~ — l~) 



.a 

'Bafi 



(27T) d 
1 



k 2 k~ 

7"6(p + q — l) 



x 



*%£Hn 



]ir 2 ia s N c C F 
1 



-p / -d+/ + m b 

d d k 5(l'~ - I- -k 



(2tt)< 



k 2 



1(1)1" 



I 



lyb(P + q-l) 



In the case of the .B-meson only the diagrams in fig. 13.9( d). (e) give rise to UV- 
poles. Those diagrams however do not lead to a mixing of Oq and the evanescent 
operators and we do not have to deal with evanescent operators. 

First let's have a look at the convolution integral which belongs to the diagram 
in fig. E£S](f) B 



(0i + o 2 )% U{) 



(AnfialN.Cl^ j 



d d k 



(27r) d 2(k + I) -pk 2 



x 



7 



(k + iy 



7 / W(l-7s) 



m b 



k 2 - 2k ■ (p + q - I) 



It 



2p~ 

y 



(3.68) 



—liflu (1 -7b)- 

yy 



In leading power ( 13. 68ft is identical to the contribution of the two diagrams shown 
in fig. 13.101 which is given by: 



(An) 2 ia 2 s N c C 2 F 



d d k 



(2n) d k 2 (k + l-xp) 2 



x 7 , no (! - 7s) ,.., 0> . — 1 : n"7r 



® 7m 



(A; + /) 2 
y4+xi 



' (yq + xp — k — I) 



k 2 -2k- (p + q-l) ' T 

I yd + xp'- ft- / 

7^ - 7^ — , ~ r — RT7m 



(yg + xp — k — iy 



2 In order not to confuse the reader I stress that the symbols 
equations are meant in terms of (|3.2|) and not of (|3.46|) . 



I-75). (3.69) 
in this and the following 
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Figure 3.10: Two diagrams which correspond in leading power to the contribution 
of the 5-meson wavefunction (13.681) . 

<? ffl . < f(xp) 



_q(xp) 



Figure 3.11: a s contributions to the form factor 

In ( 13.691) the leading power comes from the region where k is soft. In this region of 
space the integrand gets the form of the integrand in (13.681) . so both contributions 
cancel. As we did not include the diagrams of fig. 13.101 in the last section we can 
skip the contribution of ( 13.681) here. 

The remaining contributions are the diagrams in fig. 13.9( d) and (e). Together 
they read: 

(Oi + £>2>£ ),(d) ' {e) = (3.70) 
aZN&^-f&fil - 75) ® ( — 9,u ~ 4tmt0 x 

& \ y yy J 



21 n ^ ( 4 + 21nO-2f— + 21n^) +4- — -21n 2 e). 



euv rn b J \e m m b 

3.3.6 Form factor contribution 

Finally we have to calculate the contribution of (I3.42j) . It is given by 

Aormfac, = ® ® ^ ® 41 fact . ® = (3.71) 

^f {1) ' u UvqW - tbMvq) T {1 \y). 

The form factor f^' u is the a s correction of the matrix element 

(q'(xp)q(xp)\qY(l ~ l 5 )b\b(p + q - l)q'(l)), (3.72) 
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where q' e and q e are the spinors with the flavour quantum numbers of the emitted 
pion and T^\y) is given by [TUj : 



T (1) (y) 



-6['± + ln^ 



-18 + 3 



l-2y 
V 



\ny — in 



(3.73) 



2U 2 (y)-\n 2 y 



2\ny 



- (3 + 2ivr) lny - (y <-> y) 



We get /W'" by evaluating the diagrams in fig. 13.111 Using the notation of (13.21) we 
finally obtain: 



A 



1 



formfact. 



/ 



A^cCj,— 7 M ® TmT^^ 1 - 7 5 ) - 7"(1 - 7s) 



£2^+ c/+ 1 



7,(i-7fe)r (1) (i/). 



-7 M ) ® 
(3.74) 



Chapter 4 
NLO results 



4.1 Analytical results for T{ and T\ 

After the analysis of the last chapter we finally obtain the 0{a 2 s ) results for the 
hard spectator scattering kernels Tl\ which are defined by (I2.45p . Those expres- 
sions appear in convolution integrals with wave functions, where x, y and £ are 
the integration variables as defined in f)2.45p . The ultraviolet divergences are renor- 
malised in the MS-scheme. The infrared divergences drop out after subtracting the 
wave function contributions from the amplitude. The infrared finiteness together 
with the finiteness of the convolution integrals ensures that the framework of QCD- 
factorization works at this order in a s . 

The explicit 0(a 2 s ) contributions for T x n 2 read (see next page): 
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ImT« = - 2 ^4f >< (4-4) 
AN 2 m B ^ v 1 

x 2 2x 1 1 



+ 7 =vJ + + — lny 



(y — x) 3 — x) 2 y — x yx / 

x 2 2x 1 x 



+ 



— x) 3 (y — x) 2 y — x (y — x)x yx 
xlny x 3 3 



lnx 



(y — x)x — x) 2 2(y — x) 2yx 



The a 2 corrections of the hard spectator interactions have already been calcu- 
lated in [T3J Ej- However both of these calculations have been performed in the 
framework of SCET, while my result is a pure QCD calculation. In order to com- 
pare fl4.1l) - fl4.4l) to [131 EH] we have to take into account the definition of A#. The 
SCET calculation naturally uses the A#, defined by the HQET field for the 6-meson, 
while I define Xb by QCD-fields. Those two definitions differ at 0(a s ), which has 
been discussed in appendix O The difference in the logarithmic moments of the 
B-meson wave function does not play a role, because these moments occur first at 
NLO. Using (IC.14p I figured out with the help of a computer algebra system, that 
(HHJ-dUlD reproduce the results of [HI Ej. 

The main difference between the present QCD calculation and the framework 
of SCET is the way how to make the expansion in While in the QCD 

calculation this expansion takes place at the level of the amplitude and of Feynman 
integrals, in SCET the Lagrangian is expanded in powers of AqcD/ m b- This leads to 
the fact that the structure of the calculation of [TBI [T4"] is completely different from 
the present calculation such that comparing intermediate results like single Feynman 
diagrams is not possible. So it is allowed to state that the analytical coincidence of 
the present result with the results published before gives more independent test of 
[T3~l 114"] than than a SCET calculation could provide. 



4.2 Scale dependence 

It is instructive to have a closer look to the scale dependence of f l4.ip -( r4~4l) . As stated 
in section 12.2.41 the scale dependence of A 11 vanishes. In our case we can prove this 
up to 0(a 2 ) i.e. 

ik* = 0{a ' 1 (45) 

We need the scale dependence of the following quantities, which I took from [T5] : 

' -d = (12C N C 2 + (12C F - 6C G )Ci) (4.6) 



dhx/i in 
(i ty 

—— C 2 = -± (UCjfd + (12C F - 6C G )C 2 ) (4.7) 
a In jj, Air 

ik a - = (4 ' 8) 
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where (3q = llC ° ^ nf ° N and rif = 5 the number of active flavours. Regarding the 
wave functions we need the scale dependence of their convolution integrals with the 
LO kernels: 



d 



dln/j J Q x 
d r 1 d£ 



dx a s f 3 + 21nx 

—(p n {x,n) = —C F / dx — q> v [x,n) 

it J 2x 

-0m(e,/i) = y-C f I ^(41n£ + 6)0 m (£^). 



(4.9) 
(4.10) 



(14. 9p can be obtained using the renormalisation group equation (RGE) for light- 
cone distribution amplitudes, which can be found in [HI El]. We get (14.101) from |35j . 
where the RGE for the S-meson light-cone wave function, defined in the framework 
of HQET, is given. We get the RGE for the pure QCD defined wave function by 
matching the nonlocal heavy to light current with QCD. The matching coefficient 
is given in Appendix O 

In the case of B° — > 7r + 7r~ we obtain for the /x-dependent part of A 11 : 



A n (B° - 



7T + 7T~ 



Co 



a\C F 



G f 1 1 

i^KflfB J dxdyd^(j) 7T (x)(f) 7T (y)(f)Bi(0—^ x 

fln/i^l2CjvCi 



C 2 (y C n + C F (30 + 4 In x + 4 In y + 4 In - y G G 



(4-11) 



where the ellipsis (...) stands for ^-independent terms. Using f l4.6p - fl4.10p it is easily 
seen that (14.51) is fulfilled. In the case of B° — > tt 7t° one just has to interchange C\ 
and C 2 . 



4.3 Convolution integrals and factorizability 



By looking at the hard scattering kernels of (14. ip - (14.41) it is not obvious that there 
remain no singularities in the convolution integrals over wave functions (12.451) . It is 
however possible to perform the integration analytically, which proves the factoriz- 
abilty. 

Regarding the 5-meson wave function we will obtain the result in terms of the 
quantities A# and X n , which are defined in (I2.18P and (I2.19p . The 7r-meson wave 
function is given in terms of Gegenbauer polynomials: 



4>Tr(x) = 6XX 



1 



»=1 



(4.12) 



Because of the symmetry properties of the pion we set a^n-i = 0? furthermore we 
neglect a^ n for n > 2. So we need only the second Gegenbauer polynomial which is 
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given by: 



Cf /2) (x) 



(4.13) 



Using (E35D we get for the NLO of A spcct : 
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2 i fnfB^ m B 



C N 120 In — -48Ai + 152 



» , l *_2nx 1566 1008 ^ 0^7 2 
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-338tt 
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l35~ 
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+ i -292tt 



1482 
~35 



-7T 



(4.14) 



and 



l (2) 



A * ' = ry 

spect. 2 4^7-2 



, u 1467 252 w , 9 
108 In + — — + — C 3 - 6tt 2 + it 54vr 
m b 10 5 

H 40281 2926£ 



12 



-7T 



+ a£ 216 In — + 



m b 140 



35 



<(3) - 112tt 2 + i 1187T 



108 
135" 



-7T 



(4.15) 



The finiteness of the above equations proves factorization of the hard spectator 
interactions at NLO. 
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CKM-parameters 

Kd m Kd Kb i I Kb/Kb I 7 

0.974 -0.23 0.041 0.09±0.025 (70±20)deg 

Parameters of the £>-meson 

m B f B [36] Jp^ W\ Ax m r B ± r B o 

5.28GeV (210±19)MeV 1.56 ±0.17 -3.2 ± 1 11 ±4 1.67 ps 1.54 ps 



a 2,ll 



Parameters of the 7T-meson 
/£ * PSl El SO] /„. m n a\ a\ [HI S2] 

0.28 ±0.05 131MeV 130MeV 0.3 ± 0.15 

Quark and W-boson masses 
m b (m b ) rn c (m b ) m t (m t ) [TD] M w 

4.2 GeV (1.3±0.2) GeV 167 GeV 80.4 GeV 

Coupling constants 
A (5) r 

A MS Uf 

225 MeV 1.16639 x 10~ 5 GeV~ 2 

Table 4.1: Input parameters, which were used in the numerical analysis. All param- 
eters given without explicit citation can be found in [43]. Unless otherwise stated 
scale dependent quantities are given at fi — lGeV. 

Including the contributions of A y sp ' ect 1 and A y sp ' cct 2 the quantities a^n and 
defined in f l2~3B and (g^gj) are 

- 1 tr j.r a® j.r a® \ 

"1,11 — f fBn 2 _ l ( ^2^1 sp ect. 1 ' < - / 2^1 sp ect. 1 ^l^spect. 2 J 
JttJ+ Wis 

a 2, n = f fBn 2"(^l^spect. 1 ± Cl^spect. 1 ± ^^spect. 2): (4.16) 

JttJ+ m B 

A<£„. 1 = ~'^ TO ' M * mB 9(l + a' 2 f. (4.17) 

4.4 Numerical analysis 
4.4.1 Input parameters 

For my numerical analysis I use the parameters given in table 14.11 The decay 
constant f B and the ratio ^ x have been obtained by QCD sum rules in [SB] and 

[37] respectively. The logarithmic moments Ai and A2 where calculated in [13] using 
model light-cone wave functions for the S-meson [HJ H5j HfJ H7J. For the form 
factor f^ 71 I use the value from [40], which has been obtained by QCD sum rules. 



where 



4.4 Numerical analysis 
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This value is consistent with quenched and recent unquenched lattice calculations 
[38. 39J. The first Gcgenbaucr moment of the pion wave function is zero due to 
G-parity while the second moment has been obtained by lattice simulations jHJH2]. 



4.4.2 Power suppressed contributions 

In our numerical analysis we include two contributions which are suppressed in lead- 
ing power but numerically enhanced: The twist-3 contributions and the annihilation 
topologies. 

The twist-3 contributions come from the twist-3 contributions of the pion wave 
functions and are suppressed by the factor 

, , 2ml 



m b (p)(m u (p) + m d (p)) 



(see (12.331) ). which is formally of 0{KQco/mb) but numerically about 0.9. These 
contributions modify the hard spectator scattering function H n7T ( I2.50P to [TU] 

H " = Mr* / f*»<« t (Jo ' d i My) + M (4 - 18) 

where Xjj is parametrised by: 

X H ={l+p H e l ^)\^ (4.19) 

with 

< p H < 1. (4.20) 

For numerical calculations we set 

In ^-2.4. (4.21) 

Aft 

The annihilation contributions are parametrised in (12.341) - (12.371) . Analogously 
to (I4.19P Xa is parametrised by 

X A = (1 + p A e»*) In ^ (4.22) 

with 

< p A < 1. (4.23) 



4.4.3 Amplitudes a\ and a,2 

The amplitudes a\ and a 2 are defined in (12.311) . Their hard scattering parts a^u and 
a 2j n, i.e. the parts of a\ and a 2 , which contribute to A 11 (see (12.291) ). are plotted in 
fig. 14.11 as functions of the renormalisation scale \x. The strong dependence on p of 
the real part of LO is reduced at NLO. Taking the twist-3 contributions into account 
does not increase the /i-dependence too much. The imaginary part, which occurs 
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Figure 4.1: Contribution of the hard spectator corrections to a\ and a,2 as a function 
of the renormalisation scale /i. The upper two figures show the real part, where the 
LO is given by the dashed line, while the sum of LO and NLO is shown by the thick 
solid line. The twist-3 corrections are included in the graph given by the thin solid 
line. The third figure shows the imaginary part, which occurs first at 0{a 2 s ). So no 
distinction between LO and NLO is made. 
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first at NLO, is strongly dependent on the renormalisation scale. An appropriate 
choice for the scale of the hard scattering amplitude is the hard collinear scale 

/i hc = 1.5GeV (4.24) 

In the following numerical calculations we will evaluate a^n and a 2j ii at fih c . The 
vertex corrections A 1 will be evaluated at 

fi b = 4.8GeV. (4.25) 
Using the parameters of table 14.11 we obtain 

ai =1.015 + [0.039 + 0.018z] v + [-0.012] tw3 + [-0.029] LO 

+ [-0.010 - 0.03H] NLO 
a 2 =0.184 + [-0.171 - 0.080z] K + [0.038] tw3 + [0.096] LO 

+ [0.021 + 0.045i] NL o- (4.26) 

These equations are given in a form similar to (61) and (62) in [13J. The first 
number gives the tree contribution, the vertex corrections are indicated by the label 
V, the twist-3 contributions, which come from the last part of (14.181) . are labelled by 
tw3. The hard scattering part is separated into LO and NLO. The hadronic input 
parameters I used are slightly different from [13] and in contrast to [13] I evaluated 
all quantities, which belong to the hard scattering amplitude, at the hard collinear 
scale /Xhc- This is why the values I get for a\ and a 2 are different from [T3~| . 

The hard scattering amplitudes a^n and 02,11 together with their numerical errors 
read: 

a 1;II = - 0.051 ± 0.011(param.)l^(scale) ± 0.012(tw3) 

+ [-0.031 ± 0.008(param.)+°;°^ ( sca ie) ± 0.012(tw3)]i 

a 2 ,n =0.15 ± 0.03(param.)l°;^(scale) ± 0.04(tw3) 

+ [0.045 ± 0.012(param.)lg;g^(scale) ± 0.038(tw3)]z. (4.27) 

The first error comes from the error of the input parameters in table 14. 11 The scale 
uncertainty is obtained by varying fi^ c between lGeV and 6GeV. The error labelled 
by tw3 gives the error of the twist-3 contribution, which is obtained by varying pn 
between and 1 and 4>h between and In. Within the scale uncertainty (14.271) is 
compatible with [T3] . 

It is important to remark that the result I obtained in QCD comes with formally 
large logarithms lnAqcD/^fe- in contrast to the SCET calculation of [T31 [H] it 
is not possible to resum these logarithms by a pure QCD calculation. Without 
resummation, however, these logarithms might spoil perturbation theory. Regarding 
this fact it is the more important that the error arising from the scale uncertainty 
in (14.271) is small enough for perturbation theory to be valid. 



4.4.4 Branching ratios 

The dependence of the CP-averaged branching ratios on the hard collinear scale is 
shown in fig. 14.21 It is obvious that the NLO corrections reduce this dependence 
significantly. 
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Figure 4.2: CP-averaged branching ratios as functions of the hard collinear scale /^hc 
in units of 10 -6 . In the graph with the dashed line only the leading order of the hard 
spectator scattering is contained, while in the solid line hard spectator scattering is 
taken into account up to NLO. 
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From the parameter set in table |4~T1 we obtain the following CP-averaged branch- 
ing ratios 

10 6 BR(5+ -> tt+tt ) = 6.05l?;|(had.)ll;l(CKM)l°;^ (scale) ± 0.27(sublead.) 

10 6 BR(5° -> tt+tT) = 9.41^ ; ^(had.)^(CKM)+^^(scale)tJ;^(sublead.) 

10 6 BR(5° -> 7r°7r°) = 0.39l°;^(had.)i^?(CKM)i°;^(scale)l^°(sublead.). 

(4.28) 

The origin of the errors are the uncertainties of the hadronic parameters and the 
CKM parameters, the scale dependence and the subleading power contributions, i.e. 
twist-3 and annihilation contributions. The error arising from the scale dependence 
was estimated by varying //& between 2GeV and 8GeV and //he between lGeV and 
6GeV. If we compare (14. 28ft to the experimental values [7]: 

10 6 BR(B + -> tt+tt ) = 5.5 ±0.6 
10 6 BR(S° -»■ tt+tt - ) = 5.0 ±0.4 

10 6 BR(£° -> ttV) = 1.45 ±0.29 (4.29) 

we note that BR(5+ — > tt+tt°) is in good agreement with the data. This quantity 
is almost independent of 7. The other branching ratios, which come with large 
errors, depend strongly on 7. This dependence is shown in fig. 14.31 The light-grey 
band gives the uncertainty that is defined in the same way as the errors in f)4.28p . 
where different errors are added in quadrature. The solid inner line gives the central 
value. The experimental values are represented by the horizontal band, whereas the 
vertical band gives the value of 7. It is obvious that the errors of the branching 
fractions are too large for a reasonable determination of 7. 

For B + — > tt+tt° and B° — > ti + tx~ Q CD-factorization is expected to work well, 
because at tree level Wilson coefficients occur in the so called colour allowed combi- 
nation d + C 2 /N c ~ 1, while B° -»• tt tt° comes at tree level with C 2 + Ci/N c ~ 0.2 
such that subleading power corrections are expected to be more important. On the 
other hand there are big uncertainties in the parameters occurring in the combina- 
tions \V u b\f+ w , fB^\ B an d a 2- m SB] and [T3] these parameters were fitted by the 
experimental values (f4T29|) of BR(B + -> tt+tt ) and BR(B° -> tt+tt"). Setting 



a£(lGeV) = 0.39 (4.30) 



leads to 



' default 



\Vub\f** - 0.80 flKftl/fO 

fB . n on ( fB 

\ 

SB / default 



f^ B • 289 (4 - 31) 



This leads to the following branching ratios: 

10 6 BR(£ + -> tt+tt°) = 5.5 ± 0.2(param.)+°J(scale) ± 0.6(sublead.) 
10 6 BR(5° -> tt+tt-) = 5.0+°|(param.)+ ;9(scale)+°;9(sublead.) 
10 6 BR(£° ttV) = 0.77±0.3(param.)l°;2(scale)i^(sublead.). (4.32) 
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Figure 4.3: CP-averaged branching ratios as functions of the CKM-angle 7 in units 
of 10 -6 . The light-grey band gives the uncertainty from the errors of table |4~T1 and 
from the twist-3 and the annihilation contributions. The solid inner line gives the 
central value. The horizontal dark-grey band gives the experimental value according 
to [7] and the vertical grey band gives the value of 7 from table H~T1 within the error 
ranges. 
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The uncertainties of the quantities that occurred in (14.301) and (14.311) have not been 
considered in the estimation of the errors in (I4.32p . The B° — > 7r°7r° branching ratio 
obtained in ( 14.321) is compatible with the value obtained in [13]. Though it is too 
low, due to the theoretical and experimental errors it is compatible with (14. 29ft . 

There are two different sources of errors. On the one hand for errors that are 
due to uncertainties of input parameters and the renormalisation scale there is at 
least in principle no lower limit. On the other hand errors arising from sub leading 
power corrections, i.e. twist-3 and annihilation contributions, cannot be reduced 
in the framework of QCD-factorization. Fig. 14.41 shows the branching fractions of 
B° — > 7r + 7r~ and B° — > 7t°tt as functions of 7. The errors arising from subleading 
power contributions are represented by the dashed lines inside of the light-grey error 
band. While in the case of B° — > 7t + tt~ this remaining error might be small enough 
for non-trivial phenomenological statements about 7, in the case of B° — > tt 7t there 
remains an error of about 30%. 
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Figure 4.4: CP-averaged branching ratios as functions of the CKM-angle 7 in units 
of 10 -6 with the input parameters (I4.30p and (j4.3ip . The dashed lines inside of the 
light-grey band give the error coming from subleading power contributions, while 
the dashed lines at the border of the grey bands are included to lead the eye. The 
meaning of the other curves and bands is the same as in fig. 14.31 besides the fact 
that the parameters occurring in (I4.30p and (14.311) were not included in the error 
estimation. 



Chapter 5 
Conclusions 



In the last decades B physics has proven a promising field to determine parameters 
of the flavour sector with high precision. It is expected that in the next few years 
the angles a and 7, which are directly connected to the complex phase of the CKM 
matrix, will be measured with an accuracy at the percent level. Furthermore the 
discovery of physics beyond the standard model will be possible. 

On the theoretical side QCD factorization has turned out to be an appropriate 
tool to calculate B decay modes from first principles, because it allows for systematic 
disentanglement of the perturbative physics and the non-perturbative physics. The 
present calculation showed that the hard spectator scattering amplitude factorizes 
up to and including 0(a^), i.e. all infrared divergences cancel and there are no 
remaining endpoint singularities. The former point is obvious after the explicit 
calculation of T n and the latter point was shown by evaluating the convolution 
integral (12.11) analytically. The explicit expressions for the hard spectator scattering 
kernel f l4.1l) - fl4.4p confirmed the result of p31 [TJ]. So they are also a confirmation 
that the leading power of the amplitudes can be obtained by performing the power 
expansion at the level of Feynman integrals rather than at the level of the QCD 
Lagrangian using an effective theory like SCET, which was done in [T31 [T4"] . 

The main challenges in the evaluation of Feynman integrals, which were made 
possible with the help of tools like integration by parts identities and differential 
equation techniques, were due to the fact that the Feynman integrals came with 
up to five external legs and three independent rations of scales. Many steps in the 
calculations of section 13.21 might look like cookery. However I dare say calculating 
Feynman integrals is cookery. 

One motivation to calculate the 0{a 2 s ) corrections of the hard spectator interac- 
tions separately is the fact, that the LO of this class of diagrams starts at 0(a s ) such 
that in order to fix the scale we need the NLO correction. The numerical results of 
section 14.4.31 show that the NLO reduces the scale dependence significantly. This is 
even more important with respect to large logarithms that arise because of the fact 
that next to the m^-scale also the hard-collinear scale ^/AqcD m fe enters the hard 
spectator scattering amplitude. In contrast to the effective theory ansatz the QCD 
calculation of this work does not allow the resummation of these logarithms. This is 
why it is a crucial point, that the NLO is numerically important but small enough 
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for perturbation theory to be valid. 

Next to the scale dependence a main source of uncertainty is due to the fact that 
we do not know hadronic quantities well enough. This might be improved in the 
next few years by lattice calculations and even determination of the hadronic input 
parameters in experiment. Also a better control of power corrections would allow 
to obtain much more precise predictions from QCD factorization. 

Finally it is important to note that the present calculation is not the complete 
order a 2 s result as the contributions of penguin contractions and the effective penguin 
operators where not considered in this thesis. Actually they play a dominant role 
in the branching ratios of B — > Ktt and CP asymmetries of B — > nn and should 
be taken into account in phenomenological applications. While writing down this 
thesis the order a 2 s of these contributions has been recently published in [35]. Also 
the 0(a 2 s ) corrections of T 1 were not part of this thesis. These contributions have 
been calculated in [IT], [12] . 

So the calculation of the present thesis is a small but very important tessera in 
the mosaic of theoretical i?-physics. 



Appendix A 

CAS implementation of IBP 
identities 

A.l User manual 

This section will give an introduction how to use my Mathematica packages lorentz . m 
and ibp . m. These packages use the rules of section 12.4.11 and the algorithm of [22J . 
You can download these files from 

http : / /www . theorie . physik . uni-muenchen . de/~pilipp 

I assume that these files are located on your hard disk in the directory path. After 
you have started your Mathematica notebook with the two lines 

<<path/lorentz .m; 
<<path/ibp.m; 

you have to set some variables. Because my program distinguishes between Lorentz 
vectors and scalars we have to define which variables are of the type vector. This is 
done with the function 

AddMomenta[pl , . . . ,pn] 

which defines the variables pi , . . . ,pn to be of the type vector. The function 

RemMomenta[pl , . . . ,pn] 
removes the attribute vector from pi, . . . ,pn and 

ShowMomenta[] 

gives list of all vector variables. Per default the variables p, q and 1 are defined to 
be vector variables. 

The syntax of defining scalar products is the same as in Tracer [50] • The OnShell- 
command 



0nShell[on,{pl,0},{p2,p3,m}, . . .] 
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defines the scalar products pi ■ pi = and p2 ■ p3 = m. By default there are the 
following definitions: 

OnShell [on , {p , 0} , {q , 0} , {1 , 0} , {p , q , 1/2} , {p , 1 , xi/2} , {q , 1 , theta/2}] 

To undo the onshell definition use the flag off instead of on. 

An integral of the form (12.511) contains the set of momenta {pi, . . . ,p n }, which 
are in general linear combinations of basis momenta e.g. {0,p M ,p /1 + yq' 1 } where 
the basis momenta are {p^,^}. To tell Mathematica which variables are the basis 
momenta we have to define the variable MomBasis. In our example we set: 

MomBasis = {p, q}; 

After this definition the function 



ExternalMomenta[pi, . . . ,p n }; 



has to be called to tell Mathematica that p\, . . . ,p n are the momenta which appear 
in the Feynman integrals. In the above example: 

ExternalMomenta[0, p, p + y * q]; 

Feynman integrals are represented by the function FInt. This function will be 
simplified applying rule [T], [2] and rule [3] of section [2.4. II After the call of the function 
ExternalMomenta[pi, . . . ,p n ] the momentum is represented by the position % at 
which it appears in the argument list. So the integral ( 12. 5 1 f> is represented by 

FInt[{{n, Ml mi}, . . . , {i t , M 2 tl m t }}, {{h, M\, mi}, . . . , {i tt , M u 2 , m u }}, 

{{jl,8i},...,{jl,Si}}] 

Because most integrals do not have propagators of the form k ■ p + M 2 , the second 
argument of FInt can be dropped such that 

FInt [{fa, M 2 , my}, . . . , {i t , M t 2 , m t }}, {{j u Sl }, . . . , {j h s t }}} 

represents the integral 

d d k s™ 1 . . . s[" 
(27r) d D™ 1 . . . D 

For example: We want to represent the integral 



mt 
t 



d d k k • (p + yq) 

(2n) d k 2 (k + p) 2 (k+p + yq) 2 ' 



(A.1) 



After the above call of ExternalMomenta[0,p,p+y*q] the momenta 0, p, p+y*q 
are represented by the numbers 1, 2 and 3 respectively. So ( 1A.1I) is represented by 



FInt[{{l, 0, 1}, {2, 0, 1}, {3, 0, 1}}, {{3, 1}}] 



A.l User manual 
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which is transformed into 

-1/2 * y * FInt[{{l, 0, 1}, {2, 0, 1}, {3, 0, 1}}, {}] 

because rule [D and rule [3] are applied and scaleless integrals vanish in dimensional 
regularisation. 

The identities from rule|4]to rule[6]are created by the function IBP. This function 
takes three or five arguments and is called by 

IBP[Denoml, Denom2, lal, la2, lb] 

or 

IBP[Denom, la, lb] 

which is a shortcut for 

IBP[Denom, {}, la, {}, lb]. 

The first two arguments of IBP Denoml and Denom2 are lists that take the form 
{{ii, M x 2 }, . . . , {i t , M 2 }} and {{i, Mf}, . . . , {i u , M 2 }} respectively and describe the 
topology i.e. they tell Mathematica to create identities of integrals which contain the 
propagators Di, . . . , D t and D±, . . . , D u respectively. The powers m; and fhi to which 
the propagators Di and Di have to appear, are given by the next two arguments lal 
and la2. These are lists whose elements are of the form {n, m}: For all integrals 
with n different propagators of the form D and D respectively identities of the form 
(12.591) . (I2.60p and (12. 62ft are created for all integrals, where Y^k=\{ m k — 1) < m and 
Ylk=i(™k — I) < m respectively. The third argument lb has the same form as lal 
and la2 and tells the program how many scalar products of the form k ■ p should 
be in the numerator. Here n stands for the number of different propagators of the 
form D and D and the integrands have to fulfil the condition X^l=i n jk — m - F° r 
all of the lists lal, la2 and lb the default value for m is e.g. {{2, 1} , {3,0}} and 
{{2,1}} lead to the same result. 

In our example a convenient call of IBP would be 

subslist = IBP[{{1, 0}, {2, 0}, {3, 0}}, {{2, 1}, {3, 0}}, {{2, 1}, {3, 0}}]; 

which is equivalent 

subslist = IBP[{{1, 0}, {2, 0}, {3, 0}}, {{2, 1}}, {{2, 1}}]; 

because the default value of the powers of propagators is 0. The command 

FInt[{{l, 0, 1}, {2, 0, 1}, {3, 0, 1}}, {{3, l}}]//.subslist 

reduces our integral to 

-(1 - 2 * e) * FInt[{{l, 0, 1}, {3, 0, 1}}, {}]/2/e 

where e = (4 - d)/2. 



66 



A. CAS implementation of IBP identities 



A. 2 Implementation 

In this section I will describe in detail how to implement the rules of section l2.4.1l in 
Mathematica . I will follow the algorithm described in [22], the reader is expected 
to be familiar with this algorithm. 

First of all we have to tell our computer algebra system how to handle Lorentz 
vectors. We will write all the definitions into the file lorentz. m. It proves to be 
useful to distinguish between Lorentz vectors and scalar variables. In a CAS which 
does not know about type declarations of variables this is done by putting all the 
vector variables in a list we call MomList. So the first part of the file lorentz. m 
looks like this: 

BeginPackage [ "LORENTZ"' ]; 
OPattern variables*) 

Unprotect [a,b,c,mom,moml,mom2,ip,rest] ; 
Clear [a , b , c , mom , moml , mom2 , ip , rest] ; 
Protect [a,b,c,mom,moml,mom2,ip,rest] ; 

Unprotect [d,e] ; 
Clear [d,e] ; 

d = 4 - 2*e ; (*Dimension*) 
Protect [d,e] ; 

Unprotect [MomList] ; MomList := O; 

(*List of variables which are defined to be momenta; 
all other variables are handled as scalars*) 
Protect [MomList] ; 

Unprotect [AddMomenta] ; 
Clear [AddMomenta] ; 

AddMoment a [ql i st ] : = 

(Unprotect [MomList] ; 

MomList = Union [MomList , {qlist}-]; 

Protect [MomList] ; ) 
Protect [AddMomenta] ; 

Unprotect [RemMomenta] ; 
Clear [RemMomenta] ; 

RemMoment a [ql i st ] : = 

(Unprotect [MomList] ; 
MomList = Complement [MomList , -[qlist}-]; 
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Protect [MomList] ; ) 
Protect [RemMomenta] ; 

Unprotect [ClearMomenta] ; 
Clear [ClearMomenta] ; 
ClearMomenta [] : = 

(Unprotect [MomList] ; MomList = O; Protect [MomList] ; ) 
Protect [ClearMomenta] ; 

Unprotect [ShowMomenta] ; 
Clear [ShowMomenta] ; 
ShowMomenta [] := Return [MomList] ; 
Protect [ShowMomenta] ; 

After we have defined which variables are used as pattern variables and we 
have set the dimension variable d = 4-2*e, we introduce the protected list vari- 
able MomList which is manipulated by the functions AddMomenta, RemMomenta, 
ClearMomenta and ShowMomenta. 

In the next step we define the function IsVector which tells us if a variable is 
of the type vector (i.e. it is contained in MomList) or scalar: 

Unprotect [IsVector] ; 
Clear [IsVector] ; 

IsVector [a_ + b_] := IsVector [a] I I IsVector [b] ; 
IsVector [a_ b_] := IsVector [a] I I IsVector [b] ; 

IsVector [a_ /; MemberQ [MomList , a]] := True; (*that's the point*) 
IsVector [a_] := False; 
Protect [IsVector] ; 

Following the conventions of "Tracer" [50] we define a scalar product between 
Lorentz vectors which gets the attribute Orderless. The scalar product is defined to 
be linear, which makes the distinction between scalar and vector variables necessary. 

Unprotect [SP] ; 
Clear [SP] ; 

SetAttributes [SP, Orderless]; (*scalarproduct of two Lorentzvectors*) ; 
SP[(a_ /; ! IsVector [a] )*b_, c_] := a*SP[b, c] ; 
SP[a_ + b_, c_] := SP[a, c] + SP[b, c] ; 
SP[0, J := 0; 
Protect [SP] ; 

The function OnShell is defined as in [50] i-e. we define the scalar products 
p ■ q = ~ and p 2 = by the command OnShell [on, {p ,0} , {p,q, 1/2}] . 

Off [General : : spelll] ; 
Unprotect [OnShell] ; 
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Clear [OnShell] ; 

OnShell[flag_, list ] := 

(*Defined as in Tracer*) 
Module [{1, i}, 
1 = -[list}; 
Unprotect [SP] ; 
Switch [flag, 

on, For[i =1, i <= Length [1] , i++, 
Switch [Length [1 [ [i] ] ] , 

2, SP[l[[i]] [[1]] , l[[i]][[l]]] 

3, SP[l[[i]] [[1]] , l[[i]][[2]]] 
] 

], 

off, For[i =1, i <= Length [1] , i++, 
Switch [Length [1 [ [i] ] ] , 

2, SP[l[[i]] [[1]] , l[[i]][[l]]] 

3, SP[l[[i]] [[1]] , l[[i]][[2]]] 
] 

] 

]; 

Protect [SP] ; 

]; 

Protect [OnShell] ; 

The last function we introduce is Project. This function, applied to a linear 
combination of vector variables and a vector variable p, gives the coefficient of p 
in that linear combination. E.g. assume p and q are vector variables and x and y 
are scalars then Project [x*p+y*q,p] is simplified to x. This is the definition of 
Project: 

Unprotect [Project] ; 
Clear [Project] ; 
Project [a_+b_ ,mom_] / ; 

MemberQ [ShowMomenta[] ,mom] :=Project [a, mom] +Project [b,mom] ; 
Project [a_*moml_ ,mom2_] / ; 

MemberQ [ShowMoment a [] ,mom2] &&IsVector [moml] && ! IsVector [a] := 
a*Project [moml,mom2] ; 
Project [0,mom_] :=0; 

Project [mom_,mom_] /; MemberQ [ShowMoment a [] ,mom] :=1; 
Project [moml_ ,mom2_] / ; 

MemberQ [ShowMoment a [] ,mom2]&&MemberQ[ShowMomenta[] ,moml] :=0; 
Protect [Project] ; 

The end of the file lorentz.m is special for the calculation in this thesis i.e. we 
introduce the Lorentz vectors p, q and 1 whose scalar products fulfil our kinematical 
conditions: 



= l[[i]] [[2]] , 
= l[[i]] [[3]]; 
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(* 



other definitions 



*) 



AddMomenta[p, q, 1] ; 



OnShell[on, {p,p,0}, {q,q,0}, {p,q,l/2}, {1,0}, {p , 1 , xi/2} , 
{q,l,theta/2}] ; 



(* 



end other definitions 



*) 



EndPackage [] 

We will write the definitions which handle the reduction of Feynman integrals 
into the file ibp.m. The first part of this file covers the global variables which 
will be described in more detail when they will be used. The variable MomBasis 
which has to be defined by the user contains all the basis momenta i.e. all external 
momenta which appear in the integrals have to consist of a linear combination of 
the components of MomBasis. 

OPattems*) 

Unprotect [Denom,Denoml , Denom2 , Num , arg , 

exprl , expr2 , ip, ilp, i2p, i3p, i4p, intp,nlp ,n2p,Mp,Mlp,M2p , 
a,b,c,dp,al,bl,cl,dl,el,a2,b2,c2,d2,Ml,M2,ml,m2,mp,pl,p2,L, sp] ; 

Clear [Denom , Denoml , Denom2 , Num , arg , 

exprl , expr2 , ip, ilp, i2p, i3p, i4p, intp,nlp ,n2p,Mp,Mlp,M2p , 

a , b , c , dp , al , bl , c 1 , dl , e 1 , a2 , b2 , c2 , d2 , Ml , M2 , ml , m2 , mp , pi , p2 , L , sp] ; 

Protect [Denom , Denoml , Denom2 , Num , arg , 

exprl , expr2 , ip, ilp, i2p, i3p, i4p,nlp,n2p, intp,Mp,Mlp,M2p , 

a , b , c , dp , al , bl , c 1 , dl , e 1 , a2 , b2 , c2 , d2 , Ml , M2 , mp , ml , m2 , pi , p2 , L , sp] ; 

(*all global variables*) 
MomBasis = {p,q,l}; 

Unprotect [eqlist] ; 

eqlist = 0;(*IBP1 writes into eqlist*) 
Protect [eqlist] ; 

Unprotect [Mom] ; 
Mom = {}; 
Protect [Mom] ; 

Oglobal list of external momenta. 
This list is used by all other functions. 
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This list is set by the function ExternalMomenta (see below) , 
do not edit ! *) 

Unprotect [LinlndepMom] ; 
LinlndepMom = {}; 
Protect [LinlndepMom] ; 

OGives the position of the linearly independent momenta in Mom. 
This list is set by the function ExternalMomenta (see below) , 
do not edit ! *) 

If we define a set of external momenta, which appear in our integrals, not all of 
those will be linearly independent. To find the linearly independent ones i.e. a basis 
of our external momenta we define the function FindBasis. This function applied 
on a list of momenta gives a list of the position of the linearly independent ones. 
The algorithm is as follows: Consider a set of vectors {v 1 , . . . , v n } from which we 
want to choose a minimal subset of linearly independent vectors which form a basis 
of spanjf i, . . . ,v n }. We look for a most general solution of the equation system 



If one of the c l is not necessarily i.e. V\ can be expressed by a linear combina- 
tion of the other vectors we remove V\ from the set {v i, . . . , v n } and repeat the 
procedure until (1A.2|) gets the unique solution c l = for all i. Usually the basis 
of spanj^x, . . . ,v n } is not unique. In this case we have an ambiguity, which of the 
vectors we can remove. The function below removes this vector which comes first 
in the list {v i, . . . , v n } i.e. the vectors which are more behind have a higher precen- 
dence to be a member of the basis which is returned. This is the implementation in 
Mathematica: 

Unprotect [FindBasis] ; 
Clear [FindBasis] ; 

FindBasis [{ip_}] := {}/; (Mom[[ip]] == 0); 
FindBasis [p_List] := 

(*gives the position of the linearly independend momenta in p 
reads the variable MomBasis. If this is not unique, a vector 
p[[il]] has a higher precendence than p[[i2]] iff il>i2. The 
vector p[[i]] is an integer s.t. Mom[[ p[[i]] ]] is the 
corresponding vector*) 
Catch [ 

Module [{pc, Hist, V, C, c, i, j}, 




(A.2) 



pc = Array [Function [{i}, {Mom[[ p[[i]] ]], i}] , Length [p]]; 
(*maps every momentum onto a position in the list p*) 



Hist = Range [Length [p] ] ; 
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(*result list*) 

While [True , 
V = Array [ 

Function [{i , j}, Project [pc [ [i , 1]], MomBasis [ [j] ] ] ], 
{Length [pc] , Length [MomBasis] }] ; 
C = Array [c , {Length [pc] }] ; 

Off [Solve : : svars] ; 

C = C /. (Solve[C.V == 0, C][[l]]); 

On [Solve: : svars] ; 

For[i =1, i <= Length [C] , i++, 

(*Because we start with lowest value i=l and increase i, 
we make sure that the vectors in the begining of p are 
dropped out first, i.e. p[[il]] has a higher precendence 
than p[[i2]] iff il > i2.*) 
If [! MatchQ[C[[i]] , 0] , 

Hist = Drop [Hist, Position[Ilist , pc[[i, 2]]][[1]]]; 

pc = Drop[pc, {i}] ; 

Break [] ; 

]; 

If[i >= Length [C], Throw [Hist] ] ; 

]; 

]; 

] 

]; 

Protect [FindBasis] ; 

In the next part we define some global functions. The function MakeComb takes 
the integer arguments n and k and gives a list of all subsets of { 1 , . . . , n} with 
length k. The functions GetBasis, RepMom and rulel are well commented such 
that we do not explain them. These functions are defined such that they remember 
the value they have already calculated. So there are the commands def GetBasis, 
defRepMom and defrulel which redefine those functions and which are evaluated 
in the function ExternalMomenta. The argument of ExternalMomenta is a list 
{pi, . . . ,p n } of the external momenta which appear in the integrals. By calling the 
function ExternalMomenta[pi, . . . ,p n ] the variable Mom is set to {pi, . . . ,p n } and to 
the variable LinlndepMom a list of the position of the linearly independent momenta 
in {pi, . . . ,p n } is assigned. 

Unprotect [MakeComb] ; 
Clear [MakeComb] ; 
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MakeComb [n_ , k_] : = 

(*generates all subsets of {1, . . . , n} with k elements *) 

Catch [ 

Module [{LO, LI, i, j, 1}, 

LO = Array [{#} &, n] ; (*L0 = {{1}, {n}}*) 
LI = {}; 

For[l = 1, 1 <= k - 1, 1++, 
For[j = 1, j <= Length [LO], 

For[i = Last[L0[[j]]] + 1, i <= n, 

LI = Append [LI, Append [LO [ [j] ] , i]]; 

]; 

]; 

LO = LI; 
LI = {}; 

]; 

Throw [LO] ; 

]; 

]; 

Protect [MakeComb] ; 



Unprotect [partfrac,defpartf rac] ; 

(*This function gives a list {True , {cl , . . . , cn}} s.t. 
cl+. . .+cn=0, 

cl*Mom [ [mom [ [1] ] ] ] + . . . +cn*Mom [ [mom [ [n] ] ] ] =0 , 
cl*(mom[l] -mass [1] ) + . . . +cn* (mom[n] -mass [n] )=1 
If this not possible it gives back {False, {}} 
*) 

Clear [partf rac, defpartf rac] ; 

defpartfrac:=( 

partf rac [mom_,mass_] := 

Catch [ 

Unprotect [partf rac] ; 
Module [{res}, 
res=( 

partf rac [mom , mass] = 
Catch [ 

Module [{coef f , coef f list , solvelist , eql , il , j 1 , i} , 
Off [Solve : : svars] ; 
i = Length [mom] ; 
coef f list = Array [coef f , i] ; 
eql = 

Array [Function [{ j 1} , 
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Sum[coef f [il] *Project [Mom[ [ mom[[il]] ]], MomBasis [ [j 1] ] ] , 
{il, 1, i}] == 0], Length [MomBasis] ] ; 
eql = Append [eql, Sum [coef f [il] , {il, 1, i}] == 0]; 
eql = Append [eql, Sum [coef f [il] * ( 

SP [Mom [ [mom [ [i 1] ] ] ] , Mom [ [mom [ [i 1] ] ] ] ] -mass [ [il] ] 
),{il,l,i}] == 1]; 
solvelist = Solve [eql, Array [coef f, i]]; 
If [MatchQ [solvelist , {}] , Throw [{False , {}}] ] ; 
solvelist=solvelist [ [1] ] ; 
coef f list = coef f list /. solvelist; 
coef f list = coef f list /. coef f [_] -> 0; 
On [Solve: :svars] ; 
Throw [{True , coef f list}] ; 
] OendMod*) ; 
] OendCatch*) 

); 

Protect [partf rac] ; 
Throw [res] ; 
] ; OendMod*) 
] ; OendCatch*) 

); 

Protect [partf rac, def partf rac] ; 

Unprotect [def GetBasis , GetBasis] ; 
Clear [def GetBasis, GetBasis] ; 

def GetBasis := 

( GetBasis [Denoml_List ,Denom2_List] := 
OGets Denoml and Denom2 of the form 
{{pl,Ml~2,ml}, . . . ,{pt,Mt~2,mt}} and 

{{ql,Mql"2,mql}, . . . . ,{qr,Mqr"2,mqr» resp. and finds the 
linear independent momenta {q_kl , . . . ,q_kl} of {ql,...,qr} 
and {p_il , . . . ,p_il} of {pi , . . . . ,pt} , which are completed 
to a basis s.t. as much as possible vectors of 
{q_kl , . . . . ,q_kl} are in the basis list*) 

Catch [ 
Unprotect [GetBasis] ; 
Module[{res, pl,p2,p, bll,bl2, i}, 
pl=Array[ Denoml [[#,1]]&, Length [Denoml] ]; 
p2=Array[ Denom2[[#,l]]&, Length [Denom2] ]; 
p = Join [pi ,p2] ; 
For [i=l , i<=Length [Mom] , i++ , 
If [ ! MemberQ [p , i] , p=Prepend [p , i] ] ; 

]; 
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bll = FindBasis[p] ; 

bl2 = Array [ p[[ bll[[#]] ]]&, Length [bll] ]; 
res = (GetBasis [Denoml ,Denom2] = bl2) ; 
Protect [GetBasis] ; 
Throw [res] ; 
] (*endMod*) ; 
] OendCatch*) ; 

); 

Protect [def GetBasis, GetBasis] ; 

Unprotect [RepMom,defRepMom] ; 
Clear [RepMom, def RepMom] ; 
def RepMom : = ( 
RepMom [mom_,basis_List] := 

(*RepMom takes as arguments a momentum p and a basis {pi , . . . ,pk} . 

It gives a list {cl,...,ck} s.t. p = cl*pl+. . .+ck*pk *) 
Catch [ 

Module [{coef f , coef f list , solvelist , eql , i , il , i2} , 
Unprotect [RepMom] ; 

Off [Solve : : svars] ; 
eql = Array [ 
Sum[ 

coeff [il]*Project [Mom [[ basis [[il]] ]], MomBasis [ [#] ] ], 

{il , 1 , Length [basis] } ] == 
Project [Mom [ [mom] ] , MomBasis [[#] ] ]&, 
Length [MomBasis] 

]; 

eql = eql //. 

{(exprl_ == expr2_) /; FreeQ [{exprl , expr2}, coeff [_]] :> 

MatchQ [exprl , expr2]}; 
(*An expression of the form "x==x" in a list of equations 
is not automatically 
transformed into "True": This has to be done by hand*) 
coeff list = Array [ coeff, Length [basis] ]; 
solvelist = Solve [eql, coeff list] ; 
coeff list = coeff list /. solvelist; 
coeff list = coeff list /. coef f [_] -> 0; 
On [Solve : : svars] ; 

RepMom [p, basis] = coeff list; 
Throw [ coeff list [[1]] ]; 

Protect [RepMom] ; 
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Throw [res] ; 
] (*endMod*) ; 
] OendCatch*) ; 

); 



Unprotect [def rulel ,rulel] ; 
Clear [defrulel ,rulel] ; 
def rulel := ( 
rulel [mom_List] := 

(* rulel gives for the topology {il,...,ik} a list of 
coefficients (cl,...,ck) such that cl+. . .+ck = 1 and 
cl*Mom[[il]]+. . .+ck*Mom[[ik]] = 0. The output is of the form 
{True , {cl , . . . , ck}} or {False} if this is not possible *) 
Catch [ 

Module[{res, coef f list , coef f , eql , i,j,il,jl}, 
Off [Solve : : svars] ; 

coef f list = Array [coef f, Length [mom] ] ; 
eql = 

Array [Function [{j 1} , 
Sum[coef f [il] *Project [Mom[ [ mom[[il]] ]], MomBasis [ [j 1] ] ] , 
{il, 1, Length [mom] }] == 0], Length [MomBasis] ] ; 
eql = Append[eql, Sum[coef f [il] , {il, 1, Length [mom] }] == 1]; 
coef f list = coef f list /. Solve [eql, Array [coef f, Length [mom] ]] ; 
coefflist = coefflist /. coef f [_] -> 0; 

(*Set all the coefficients which remain after solving the 
system of equations to zero, if there is no solution of 
this system of equations it follows 

coefflist = {coeff[l], coef f [2] , coef f [3]} -> {0, 0, 0}. 
As the linear equation system is inhomogeneous , coefflist 
can only be substituted to {0, 0, 0} if there is no 
solution: 

The momenta are linearly independend or the condition 
sum coeff [i] = 1 cannot be satisfied*) 
On [Solve : : svars] ; 

If [(coefflist //. {0, r } -> {r}) == {}, 

(*first case : k~2 not reducible*) 
res={False}- ; , 

(*else - part : k~2 reducible*) 
coefflist = coef f list [ [1] ] ; 
res={True , coefflist} ; 
] ; Oendif *) 

Unprotect [rulel] ; 
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rulel [mom] =res ; 

(*remember previously calculated values*) 
Protect [rulel] ; 

Throw [res] ; 

] (*endMod*) ; 
] OendCatch*) ; 

); 

Protect [def rulel , rulel] ; 

Clear [ExternalMomenta] ; 
ExternalMomenta[L ] := 

(*sets the variables Mom and LinlndepMom*) 

Module [{i , il , j , j 1 , j2,k, eql , coef f , coef f list , 
LI ,L2,L3,L4, cmom,base, solvelist}, 

Unprotect [Mom, LinlndepMom]; 
Clear [Mom, LinlndepMom] ; 
Mom = {L}; 

(*The variable Mom has to be set. Otherwise you are not allowed 
to use the function FindBasis *) 
LinlndepMom = FindBasis [Range [Length [Mom] ]] ; 
Protect [Mom, LinlndepMom] ; 

Unprotect [partf rac] ; 
Clear [partf rac] ; 
def partf rac ; 

Odefine the function partfrac*) 
Protect [partfrac] ; 

Unprotect [GetBasis] ; 
Clear [GetBasis] ; 
def GetBasis ; 
Protect [GetBasis] ; 

Unprotect [RepMom] ; 
Clear [RepMom] ; 
def RepMom; 
Protect [RepMom] ; 

Unprotect [rulel] ; 
Clear [rulel] ; 
def rulel ; 
Protect [rulel] ; 
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]; 

Protect [ExternalMomenta] ; 

Feynman integrals are represented by the function FInt. This function accesses 
the function scaleless which gives true if the topology is scaleless. In this case 
FInt gives 0. FInt uses the functions GetBasis, RepMom and rulel to get simplified 
according to rule [IHHl 

Passarino Veltman reduction **************) 
Unprotect [Scaleless] ; 
Clear [Scaleless] ; 

Scaleless [Denom_List] := 
Catch [ 

Module [{pi, pj , i, j}, 
For [i=l , i<=Length [Denom] , i++ , 
If [Denom [ [i , 2] ] ! =0 , Throw [False] ] ; 
pi=Mom[[ Denom [[i,l]] ]] - Mom[[ Denom [[1,1]] ]]; 
For [ j =2 , j <=Length [Denom] , j ++ , 
pj=Mom[[ Denom [[j]] [[1]] ]] - Mom [[ Denom [[1]] [[1]] ]]; 
If [!MatchQ[Simplify[SP[pi,pj]] ,0] , Throw [False] ;] ; 

]; 
]; 

Throw [True] ; 

]; 

]; 

Scaleless [Denoml_List , Denom2_List] := 
Catch [ 

Module [{pi , pi , p j , Mi , i , j } , 
If [Length [Denoml] >0 , 
pi = Mom[[ Denoml [[1,1]] ]]; 
For [i=l , i<=Length [Denoml] , i++ , 
If [Denoml [ [i , 2] ] ! =0 , Throw [False] ] ; 
pi=Mom[[ Denoml [[i, 1]] ]] - pi; 
For [ j =2 , j <=Length [Denoml] , j ++ , 
pj=Mom[[ Denoml [ [j ,1]] ]] - pi; 

If [IMatchQ [Simplify [SP [pi, pj]] ,0] , Throw [False] ;] ; 

]; 

For [ j =1 , j <=Length [Denom2] , j ++ , 
pj=Mom[[ Denom2[[j , 1]] ]] ; 

If [IMatchQ [Simplify [SP [pi, pj]] ,0] , Throw [False] ;] ; 

]; 
]; 

, Oelse part *)pl=0;] Oendlf *) ; 
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For [i=l , i<=Length [Denom2] , i++ , 
pi=Mom [ [ Denom2 [ [i , 1] ] ] ] ; 

If [ ! MatchQ [tres=Simplif y [Denom2 [ [i , 2] ] -SP [pi , pi] ] , 0] , Throw [False] ] ; 

For [j=l , j<=Length[Denom2] 
pj=Mom[[ Denom2[[j , 1]] ]]; 

If [! MatchQ [Simplify [SP [pi ,pj] ] ,0] , Throw [False] ;] ; 
] (*endFor [j] *) ; 
] OendFor [i] *) ; 



Throw [True] ; 

] (*endMod*) ; 
] OendCatch*) ; 

Protect [Scaleless] ; 

Unprotect [FInt] ; 
Clear [FInt] ; 

FInt [Denom_List ,Num_List] :=0/ ; Scaleless [Denom] ; 

FInt [Denoml_List ,Denom2_List ,Num_List] :=0/ ; Scaleless [Denoml ,Denom2] ; 
(*scaleless integrals vanish*) 

FInt [Denom_List,{},Num_List] := FInt [Denom, Num] ; 

FInt [Denom_List , {a , {ip_Integer, 0}, b }] := FInt[Denom, {a, b}] ; 

FInt [Denoml_List , Denom2_List, {a , {ip_Integer, 0}, b }] := 

FInt [Denoml, Denom2, {a, b}] ; 

FInt [{a , {ip.Integer, Mp_, 0}, b }, Num.List] := FInt [{a, b}, Num]; 

FInt [Denom_List , {a , {ip_Integer, Mp_, 0}, b }, Num_List] := 

FInt [Denom, {a, b}, Num]; 

FInt [{a , {ip_Integer, Mp_, 0}, b }, Denom_List, Num_List] := 

FInt [{a, b}, Denom, Num]; 

FInt [{al , {ilp_Integer ,Mp_, i2p_Integer} , 

{i lp_Integer , Mp_ , i3p_Integer} , bl } , 

Num_List] := 
FInt [{al , {i lp , Mp , i2p+i3p} , bl} , Num] ; 

FInt [{al , {ilp_Integer ,Mp_ , i2p_Integer} , 

{i lp_Integer , Mp_ , i3p_Integer} , bl } , 

Denom_List, Num_List] := 
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FInt [{al ,{ilp,Mp, i2p+i3p},bl}, Denom, Num] ; 

Flnt [Denom_List , {al , {ilp_Integer ,Mp_ , i2p_Integer} , 

{i lp_Integer , Mp_ , i3p_Integer} , bl } , 

Num_List] := 
FInt [Denom, {al,{ilp,Mp,i2p+i3p},bl}, Num]; 

FInt [Denom_List , {al , {ilp_Integer , i2p_Integer} , 

{ilp_Integer , i3p_Integer} ,bl }] := 

FInt [Denom , {al , {ilp , i2p+i3p} , bl}] ; 

FInt [Denoml_List , Denom2_List, {al , {ilp_Integer , i2p_Integer} , 

{ilp_Integer , i3p_Integer} ,bl }] : = 

FInt[Denoml, Denom2, {al,{ilp,i2p+i3p},bl}] ; 

FInt [{al , {ilp.Integer ,Mlp_ , i2p_Integer} , 

{i3p_Integer , M2p_ , i4p_Integer} , bl } , 

Num_List]/; ilp > i3p := 
FInt [{al , {i3p , M2p , i4p} , {ilp ,Mlp , i2p} , bl} , Num] ; 

FInt [Denom_List , {al , {ilp_Integer ,Mlp_ , i2p_Integer} , 

{i3p_Integer , M2p_ , i4p_Integer} ,bl } , 

Num.List]/; ilp > i3p := 
FInt[Denom, {al,{i3p,M2p,i4p},{ilp,Mlp,i2p},bl},Num] ; 

FInt [{al , {ilp_Integer ,Mlp_ , i2p_Integer} , 

{i3p_Integer , M2p_ , i4p_Integer} , bl } , 

Denom_List, Num_List]/; ilp > i3p := 
FInt[{al,{i3p,M2p,i4p},{ilp,Mlp,i2p},bl}, Denom, Num]; 

FInt [Denom_List , 

{a , {ilp_Integer , nlp_Integer} , 

{i2p_Integer , n2p_Integer} , b }] /; ilp > i2p := 

FInt[Denom, {a, {i2p, n2p}, {ilp, nip}, b}] ; 

FInt [Denoml_List , Denom2_List, 

{a , {ilp_Integer , nlp_Integer} , 

{i2p_Integer , n2p_Integer} , b }] /; ilp > i2p := 

FInt[Denoml, Denom2, {a, {i2p, n2p}, {ilp, nip}, b}] ; 
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FInt[{{pl_, Ml_, ml_}, Denom }, {a ,{0, n.Integer} ,b }] := 

(*The entry {0, n} in the numerator denotes (k~2)~n in the 
integrand where k is the integration variable*) 

FInt [{Denom} , {a, {0, n - 1}, b}] + 
(Ml - SP[ Mom [[pi]], Mom [[pi]] ])* 

FInt [{{pi, Ml, ml}, Denom}, {a, {0, n - 1}, b}] - 
2*FInt[{{pl, Ml, ml}, Denom}, {{pi, 1}, a, {0, n - 1}, b}] ; 

FInt[{{pl_, Ml_, ml_}, Denoml }, Denom2_, 

{a ,{0, n_Integer},b }] := 

(*The entry {0, n} in the numerator denotes (k~2)~n in the 
integrand where k is the integration variable*) 

FInt [{Denoml}, Denom2, {a, {0, n - 1}, b}] + 
(Ml - SP[ Mom [[pi]], Mom [[pi]] ])* 

FInt [{{pi, Ml, ml}, Denoml}, Denom2, {a, {0, n - 1}, b}] - 
2*FInt[{{pl, Ml, ml}, Denoml}, Denom2, {{pi, 1}, a, {0, n - 1}, b}] ; 



FInt [Denom_List , Num_List] := 

(*Try to expand the denomintor into partial fractions*) 
Module [{pf ,1th}, 
Sum[pf [[2,i]]* 

FInt [ReplacePart [Denom , Denom [ [i , 3] ] -1 , {i , 3}] , Num] , {i , 1 , 1th}] / ; 
(lth=Length [Denom] ; 

pf =partf rac [Array [Denom [ [# , 1] ] & , 1th] , Array [Denom [ [# , 2] ] & , 1th] ] ; 
pf [[1]]) 

]; 

(*reduction of the HQET-Propagators*) 

FInt [Denom_List , {a , {ip_Integer , Mp_ , mp_Integer} , b } , 

{al , {ip_Integer , sp_Integer} ,bl }] : = 

FInt [Denom,{a,{ip,Mp,mp-l},b},{al ,{ip,sp-l},bl}] - 
Mp*FInt [Denom,{a,{ip,Mp,mp},b},{al ,{ip,sp-l},bl}] ; 



(*reduction corresponding rulel*) 

FInt [{a , {pl_Integer, Ml_, ml_Integer}, b }, 

{c , {pl_Integer, nl_Integer}, dp }] := 

Module [{res , unchanged, coeff list , coeff, denom, num, p, M, i, j, 1, rl}, 
(*Rule I*) 
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res/; 
Catch [ 

res= 

Catch [ 

1 = Length [{a}] + 1; OPosition of pi, Ml, ml; 

i.e. denom[[l, 1]] = pi etc. with denom see below*) 

denom = -[a, {pi, Ml, ml}, b}; 
num = {c, {pi, nl}, dp}; 

p = Array [denom [ [# , 1] ] & , Length [denom] ] ; 
M = Array [denom [ [# , 2] ] & , Length [denom] ] ; 



rl=rulel [p] ; 
If[ !rl[[l]] , 

(♦first case : k~2 not reducible*) 

Throw [unchanged] ; , 

(*else - part : k~2 reducible*) 

Throw [l/2*Sum[(KroneckerDelta[l, jl ] - rl[[2,jl]])* 
(FInt [ 

ReplacePart [denom, denom [[jl, 3]] - 1, {jl, 3}], 
{c, {pi, nl - 1}, dp}] 
+ (M[[jl]] - SP[ Mom[[ p[[jl]] ]] , 

Mom[[ p[[jl]] ]] ])* 
FInt [denom, {c, {pi, nl - 1}, dp}]), 
{ j 1 , 1 , Length [p] }] ] ; 
] ; (*endif *) 
] (*endCatch*) ; 

If [MatchQ [res , unchanged] , Throw [False] ] ; 
Throw [True] ; 
] (*endCatch*) 
] (*endMod*) ; 

FInt [{a , {pl_Integer, Ml_, ml_Integer}, b } ,Denom_List , 

{c , {pl_Integer, nl_Integer}, dp }] := 

Module [{res , unchanged, coeff list , coeff, denom, num, 
p, M, i, j, 1, rl}, 
(*Rule I*) 
res/; 
Catch [ 
res= 
Catch [ 

1 = Length [{a}] + 1; (*Position of pi, Ml, ml; 

i.e. denom [[1, 1]] = pi etc. with denom see below*) 
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denom = {a, {pi, Ml, ml}, b}; 
num = {c, {pi, nl}, dp}; 

p = Array [denom [ [# , 1] ] & , Length [denom] ] ; 
M = Array [denom [ [# , 2] ] & , Length [denom] ] ; 



rl=rulel [p] ; 
If[ !rl[[l]], 

(*first case : k~2 not reducible*) 

Throw [unchanged] ; , 
(*else - part : k~2 reducible*) 

Throw [l/2*Sum[(KroneckerDelta[l, jl ] - rl[[2,jl]])* 
(FInt[ 

ReplacePart [denom, denom [[jl, 3]] - 1, {jl, 3}], 
Denom, {c, {pi, nl - 1}, dp}] 
+ (M[[jl]] - SP[ Mom[[ p[[jl]] ]] , 

Mom[[ p[[jl]] ]] ])* 
FInt [denom, Denom, {c, {pi, nl - 1}, dp}]), 
{ j 1 , 1 , Length [p] }] ] ; 
] ; Oendif *) 
] OendCatch*) ; 

If [Mat chQ [res , unchanged] , Throw [False] ] ; 
Throw [True] ; 
] OendCatch*) 
] OendMod*) ; 



FInt [Denoml_List , Denom2_List , Num_List] : = 

(*decompose the momenta of Num into a unique set of momenta 
given by Denom2, Denoml and some futher momenta which complete 
the momenta of Denom2 and Denoml to a basis*) 
Module [{res , unchanged, base, pn, rep, i,j }, 
res/; 
Catch [ 

If [Mat chQ [Num,{}] , Throw [False] ] ; 
res=Catch [ 

base = GetBasis [Denoml ,Denom2] ; 

pn = Array [Num [ [# , 1] ] & , Length [Num] ] ; 

For[i =1, i <= Length [Num], i++, 

If[MemberQ[ base, pn[[i]] ], Continue [] ]; 
(*The momentum in the numerator already appears in 
the corresponding basis*) 
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rep = RepMom[pn[[i]] ,base] ; 

Throw [Sum [ 

rep[[j]]* 

FInt [ Denoml , Denom2 , 
Prepend [ 

ReplacePart[Num, Num[[i, 2]] - 1, {i, 2}], 
{base[[j]] ,1}]] ,{j , 1, Length [base] }] ]; 

] (*endFor [i] *) ; 
Throw [unchanged] ; 
] (*endCatch*) ; 
If [MatchQ [res , unchanged] , Throw [False] ] ; 
Throw [True] ; 
] OendCatch*) 
] OendMod*) ; 



FInt [Denom_List ,Num_List] := 

Odecompose the momenta of Num into a unique set of momenta 
given by Denom and some futher momenta which complete the 
momenta of Denom to a basis*) 
Module [{res , unchanged, base, pn, rep, i,j }, 
res/; 
Catch [ 

If [MatchQ [Num,{}] , Throw [False] ] ; 
res=Catch [ 

base = GetBasis [Denom, {}] ; 

pn = Array [Num [[#, 1]]&, Length [Num]] ; 

For[i =1, i <= Length [Num], i++, 

If [MemberQ [ base, pn[[i]] ], Continue [] ]; 
(*The momentum in the numerator already appears in 
the corresponding basis*) 

rep = RepMom[pn[[i]] ,base] ; 
Throw [Sum [ 

rep[[j]]* 

FInt [ Denom, 

Prepend [ 

ReplacePart[Num, Num[[i, 2]] - 1, {i, 2}], 
{base[[j]] ,1}]] 1, Length [base] }] ]; 

] (*endFor [i] *) ; 
Throw [unchanged] ; 
] OendCatch*) ; 
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If [MatchQ [res , unchanged] , Throw [False] ] ; 
Throw [True] ; 
] OendCatch*) 
] (*endMod*) ; 

The implementation of (12.561) is given by these two rules: 

(*these two rules should be applied last*) 

FInt[{al , {pl_Integer, Ml_, ml.Integer}, bl , 

{p2_Integer, M2_, m2_Integer}, cl }, 

{dl , -[p2_Integer , n2_Integer}, el }] : = 

Catch [ 

Module [{coeff list , denom, i, j}, 

denom = {al,{pl, Ml, ml}, bl, {p2, M2, m2}, cl}; 

Throw [1/2* ( 

FInt[{al,{pl, Ml, ml}, bl, {p2, M2, m2 - 1}, cl}, 

{dl, {p2, n2 - 1}, el}] - 
FInt[{al,{pl, Ml, ml - 1}, bl, {p2, M2, m2},cl}, 

{dl, {p2, n2 - 1}, el}] 
+ (M2 - Ml + SP[ Mom [[pi ]] , Mom[[ pi ]] ] - 
SP[ Mom[[p2]] , Mom[[ p2 ]] ]) 
*FInt [denom, {dl, {p2, n2 - 1}, el}] 
+ 2*FInt [denom, {{pi, 1}, dl , {p2, n2 - 1}, el}])]; 
] (*endMod*) ; 
] OendCatch*)/; 

MemberQ[GetBasis[{al,{pl,Ml,ml},bl,{p2,M2,m2},cl},{}] ,pl] ; 

FInt[{al ,{pl_Integer, Ml_, ml_Integer}, bl , 

{p2_Integer, M2_, m2_Integer}, cl }, Denom_, 

{dl , {p2_Integer, n2_Integer}, el }] := 

Catch [ 

Module [{coeff list , denom, i, j}, 

denom = {al,{pl, Ml, ml}, bl, {p2, M2, m2}, cl}; 

Throw [1/2* ( 

FInt[{al,{pl, Ml, ml}, bl, {p2, M2, m2 - 1}, cl}, 

Denom, {dl, {p2, n2 - 1}, el}] - 
FInt[{al,{pl, Ml, ml - 1}, bl, {p2, M2, m2},cl}, 

Denom, {dl , {p2, n2 - 1}, el}] 
+ (M2 - Ml + SP[ Mom [[pi ]] , Mom[[ pi ]] ] - 
SP[ Mom[[p2]] , Mom[[ p2 ]] ]) 
*FInt [denom, Denom, {dl, {p2, n2 - 1}, el}] 
+ 2*FInt [denom, Denom, {{pi, 1}, dl, {p2, n2 - 1}, 
el}] 

) 
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]; 

] OendMod*) ; 
] (*endCatch*)/; 

MemberQ[GetBasis[{al,{pl,Ml,ml},bl,{p2,M2,m2},cl},Denom] ,pl] ; 



The next step is to implement rule HE We define the function IBP1 which 
takes as a starting point an integral as in (12.511) and is given the two arguments 
{{il,Ml~2,ml}, . . . ,{it,Mt~2,mt}} and {{jl.sl}, . . . ,{jl,sl}} like FInt. For 
the integral defined by these arguments the identities (12.591) . (I2.60p and (12.621) are 
generated and written into the variable eqlist. 

Unprotect [IBP1] ; 
Clear [IBP1] ; 

IBP1 [Denom.List, Num_List] := IBP1 [Denom, {}, Num] ; 

IBP1 [Denoml.List , Denom2_List, Num.List] := 

Module [{expr, s, t, 1, pi, p2, Ml, M2, ml, m2, r, n, 
il, i2, i, j, idl}, 

Unprotect [eqlist] ; 

pi = Array [Denoml [[#, 1]] &, Length [Denoml] ] ; 

Ml = Array [Denoml [[#, 2]] &, Length [Denoml] ] ; 

ml = Array [Denoml [[#, 3]] &, Length [Denoml] ] ; 

p2 = Array [Denom2 [ [# , 1] ] & , Length [Denom2] ] ; 

M2 = Array [Denom2 [ [# , 2] ] & , Length [Denom2] ] ; 

m2 = Array [Denom2 [ [# , 3] ] & , Length [Denom2] ] ; 

n = Array [Num [ [# , 2] ] & , Length [Num] ] ; 

s = Sum [Num [[i, 2]], {i, 1, Length [Num] }] ; 

r = Sum [Denoml [ [i , 3]], {i, 1, Length [Denoml] }] ; 

t = Length [Denoml] ; 

1 = Length [Num] ; 



(* Identity I *) 

expr = (Collect [(d + s - 2*r) *FInt [Denoml , Denom2, Num] - 
Sum [2 *ml[[i]]*( 

(Ml[[i]] - SP[Mom[[pl[[i]]]] , Mom[[pl[[i]]]]])* 
FInt [ReplacePart [Denoml , 

Denoml [[i, 3]] + 1, {i, 3}], 
Denom2 , Num] - 
FInt [ReplacePart [Denoml , 

Denoml [[i, 3]] + 1, {i, 3}], 
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Denom2, Append [Num, {pl[[i]], 1}]]), 
{i, 1, t}]- 
Sum[m2[[i]]* 

FInt [Denoml , ReplacePart [Denom2 , 

Denom2 [ [i , 3] ] +1 , {i , 3}] , 
Append [Num, {p2[[i]] ,1}]] , 
{i , 1 , Length [Denom2] }] , 
FInt [ ]] 

); 

eqlist = {expr}; 
(* Identity II *) 

For[il = 1, il <= Length [LinlndepMom] , il++, 
i = LinlndepMom [ [il] ] ; 
expr = 

(Collect [ 
Sum[ 

n[[j]]*SP[ Mom[[i]], Mom[[ Num[[j, 1]] ]] ]* 
FInt [Denoml , Denom2 , 

ReplacePart [Num, Num[[j, 2]] - 1, { j , 2}]], 

{j, 1, 1}]- 
Sum[ 

2*ml[[j]]* 

(SP[ Mom[[i]], Mom[[ pi [[j] ] ]] ]* 

FInt [ReplacePart [Denoml , Denoml [ [j , 3] ] + 1 , { j , 3}] , 

Denom2 , Num] + 
FInt [ReplacePart [Denoml , Denoml [ [j , 3] ] + 1 , { j , 3}] , 

Denom2, Append [Num, {i, 1}]]), 
{j, 1, t}]- 

Sum[m2[[j]]*SP[ Mom[[i]], Mom [ [p2 [ [j ] ] ]] ]* 
FInt [Denoml , 

ReplacePart [Denom2, Denom2[[j, 3]] + 1, { j , 3}], 
Num] , 

{ j , 1 , Length [Denom2] }] , 

FInt [ ]]) ; 

eqlist = Append [eqlist , expr]; 
] ; OendFor [il] *) 

(* Identity III *) ; 

idl = MakeComb [Length [LinlndepMom] , 2] ; 
For[i =1, i <= Length [idl] , i++, 

11 = LinlndepMom [[ idl[[i, 1]] ]]; 

12 = LinlndepMom [[ idl[[i, 2]] ]]; 
expr = (Collect [ 

Sum[2*n[[j]]*(SP[Mom[[Num[[j, 1]]]], Mom[[i2]]]* 
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FInt [Denoml , Denom2 , 
Append [ 
ReplacePart [Num, 

Num[[j, 2]] - 1, {j, 2}], 
{il, 1}]] - 
SP [Mom [ [Num [ [ j , 1]]]], Mom[[il]]]* 
FInt [Denoml , Denom2 , 
Append [ 
ReplacePart [Num, 

Num[[j, 2]] - 1, {j, 2}], 
{12, 1}]]), 
{j, 1, Length [Num] }]- 

Sum[4*ml[[j]]*( 

SP [Mom [ [Denoml [ [j , 1]]]], Mom[[i2]]]* 
FInt [ReplacePart [Denoml , 

Denoml [[j, 3]] + 1, { j , 3}] , 
Denom2, Append [Num, {il, 1}]] - 

SP [Mom [ [Denoml [ [j , 1]]]], Mom[[il]]]* 
FInt [ReplacePart [Denoml , 

Denoml [[j, 3]] + 1, { j , 3}] , 
Denom2, Append [Num, {i2, 1}]]), 
{ j , 1, Length [Denoml] }] - 

Sum[2*m2[[j]]*( 

SP [Mom [ [Denom2 [ [ j , 1]]]], Mom[[i2]]]* 
FInt [Denoml , 

ReplacePart [Denom2 , 

Denom2 [ [j , 3] ] + 1 , { j , 3}] , 
Append [Num, {il, 1}]] - 
SP [Mom [ [Denom2 [ [ j , 1]]]], Mom[[il]]]* 
FInt [Denoml , 

ReplacePart [Denom2 , 

Denom2 [ [j , 3] ] + 1 , { j , 3}] , 
Append [Num, {i2, 1}]]), 
{ j , 1 , Length [Denom2] }] , 
FInt [ ]] 

); 

eqlist = Append [eqlist , expr] ; 
] (*endFor[i]*) ; 



Protect [eqlist] ; 
] OendModule*) ; 
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Protect [IBP1] ; 

The function IBP follows the algorithm of [22] to generate a list of the IBP 
identities. After generation a new set of identities by calling IBP1 it replaces more 
complex integrals by less complex ones. The complexity of integrals is defined in 

Unprotect [Hf 1] ; 
Clear [Hf 1] ; 

Hfl[l, Mp_] := {{Mp}}; 
Hfl[l_, 0] := {Array[0 &, 1]}; 
Hfl[l_, Mp_] := 

Ogenerate all 1 - tuples {nl , nl} s .t . nl + . . . + nl = Mp *) 

Catch [ 

Module [-[LI, L2, i, j>, 
L2 = {}; 

For[i =0, i <= Mp, i++, 
LI = Hfl[l - 1, Mp - i] ; 
For[j = 1, j <= Length [LI] , 

L2 = Append [L2, Join[{i}, Ll[[j]]]]; 

]; 

]; 

Throw [L2] ; 
] (*endMod*) ; 
] OendCatch*) ; 
Protect [Hf 1] ; 



Unprotect [BT] ; 
Clear [BT] ; 

BT [Ll_List , L2_List] := 

Ogive an ordering to lists of integers*) 

Catch [ 

Module [{i}, 

For[i =1, i <= Min [Length [LI] , Length [L2]], i++, 
If[Ll[[i]] >L2[[i]], Throw [True] ] ; 
If[Ll[[i]] <L2[[i]], Throw [False] ] ; 

]; 

If [Length [LI] > Length [L2], Throw [True] ] ; 
Throw [False] ; 

]; 
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]; 

Protect [BT] ; 



Unprotect [Verbose] ; 
Clear [Verbose] ; 
Unprotect [verbosef lag] ; 
verbosef lag=False ; 
Verbose [flag_] := 
( 

Unprotect [verbosef lag] ; 
If [MatchQ [f lag, on] , verbosef lag=True ; ] ; 
If [MatchQ [f lag, of f ] , verbosef lag=False ;] ; 
Protect [verbosef lag] 

); 

Protect [verbosef lag] ; 

Protect [Verbose] ; 

Unprotect [IBP] ; 
Clear [IBP] ; 

IBP [Denom_List , la_List, lb_List] := 
IBP[Denom, {}, la, {}, lb]; 

IBP[Denoml_List, Denom2_List, lal.List, la2_List, lb_List] := 
Catch [ 

Module [-[dummy , il, i2, j, k, 1, ml, m2, nl , n2, 
Mp, Mp2, Mdl, Md2, si, subslist, 
ell, cl2, hi, delistl, delist2, 
momlist, num, denoml, denom2, maxf, 
subsrule, al , a2, b}, 

For[i =0, i <= Length [Denoml] , i++, 

al[i] = 0; 
] ; (*Def ault value for al*) 

For[i =0, i <= Length [Denoml] +Length [Denom2] , i++, 

b[i] = 0; 
] ; (*Def ault value for b*) 

For[i =0, i <= Length [Denom2] , i++, 
a2[i] = -1; 

]; 
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a2[0] = 0; 

(*Default value for a2*) 



For[i =1, i <= Length [lal] , i++, 
al[lal[[i, 1]]] = lal[[i, 2] ] ; 

]; 

For[i =1, i <= Length [la2] , 
a2[la2[[i, 1]]] = la2[[i, 2] ] ; 

]; 

For[i =1, i <= Length [lb] , i++, 
b[lb[[i, 1]]] = lb[[i, 2]] ; 

]; 

subslist = {}; (*FInt[ . . .] -> ...*); 
For[nl = 0, nl <= Length [Denoml] , nl++, 
ell = MakeComb [Length [Denoml] , nl] ; 

For[n2 =0, n2 <= Length [Denom2] , n2++, 
cl2 = MakeComb [Length [Denom2] , n2] ; 

For[il = 1, il <= Length [ell], il++, 
For[i2 =1, i2 <= Length [cl2] , i2++, 

For[Mp = 0, Mp <= b[nl+n2], Mp++, 
momlist = Hf 1 [Length [Mom] , Mp] ; 

For[l =1, 1 <= Length [momlist] , 1++, 

num = Array [{#, momlist [[1,#]]} &, Length [momlist [ [1] ] ] ]; 
num = num//.{rl ,{il_Integer,0},r2 } -> {rl,r2}; 

For[Mdl =0, Mdl <= al [nl] , Mdl++, 
delistl = Hf l[nl, Mdl] ; 
For [Md2 = , Md2 <= a2 [n2] , Md2++ , 
delist2 = Hf l[n2, Md2] ; 
For[ml = 1, ml <= Length [delistl] , ml++, 
denoml = 

Array [Join [Denoml [[ cll[[il,#]] ]], 

{l+delistl[[ml, #]]}] &, nl] ; 
For[m2 = 1, m2 <= Length [delist2] , m2++, 
denom2 = 

Array [Join [Denom2 [ [ cl2 [ [i2 , #] ] ] ] , 
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{l+delist2[[m2, #]]}] &, n2] ; 



(*Step 8*) 
If [(MatchQ [denom2,{}]&& 

MatchQ [FInt [denoml , num] / . FInt->FInt in, 
FIntin[denoml ,num] ] ) I I 
MatchQ [ FInt [denoml , denom2, num] / . FInt->FIntin, 
FInt in [denoml, denom2, num] ], 
If [verbosef lag , Print [denoml , denom2 , num] ] ; 

IBP1 [denoml , denom2 , num] , 
Continue [] ] ; 

OCreate IBP-identities from topologies, which cannot 
be reduced by passarino veltman. *) 

Unprotect [eqlist] ; 
OStep 9(a)*) 

For[j = 1, j <= Length [eqlist] , 
eqlist [[j]] = eqlist [ [j] ]//. subslist ; 
eqlist [[j]] = Collect [eqlist [ [j] ] , 

HoldPattern[FInt[ ]] , 

Expand [ Together [#] ] &] ; 
(*substitude all the known identities into the new 
IBP identities *) 



hi = {}; 

eqlist [[j]] //. {FInt[expr ] :> ( 

hi = Append [hi , FInt [expr] ] ; 
dummy)}-; 

(*extract all the Feynmanintegrals from eqlist [[j]] 

and write them into hi *) 
If [MatchQ [hi , {}] , Continue [] ] ; 



OStep 9(b) *) 



hi = hi /. {FInt [Del. , De2_, Nu_] :> 
FInt in [Del, De2, Nu, 

Join [{Length [Del] +Length [De2] } , 

{Sum[ Del[[k,3]], {k, 1 , Length [Del] }] + 
Sum [ De2 [ [k , 3] ] , {k , 1 , Length [De2] }] } , 
{Sum [ De2 [ [k , 3] ] , {k , 1 , Length [De2] }] } , 
{Sum[ Nu[[k,2]], {k, 1, Length [Nu]}]}, 
Array [Del [[#,1]]&, Length [Del]] , 
Array [Del [[#, 3] ]&, Length [Del]] , 
Array [De2 [ [# , 1] ] & , Length [De2] ] , 
Array [De2 [ [# , 3] ] & , Length [De2] ] , 
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Array [Nu [ [# , 2] ] & , Length [Nu] ] ] ] } ; 

hi = hi /. {FInt[De_, Nu_] :> 
FIntin [De, Nu, 

Join [{Length [De] } , 

{Sum[ De[[k,3]], {k, 1, Length [De]}]}, 
{0}, 

{Sum[ Nu[[k,2]], {k, 1, Length [Nu]}]}, 
Array [De [ [# , 1] ] & , Length [De] ] , 
Array [De [ [# , 3] ] & , Length [De] ] , 
Array [Nu [ [# , 2] ] & , Length [Nu] ] ] ] } ; 

(*Give to the Feynman integrals a specific weight 
which measures the complexity of the integral*) 

maxf = 1 ; (*hl [ [maxf ] ] will be the Feynmanintegral 

with the highest complexity*) 
For[ k=2, k<=Length [hi] , k++, 
If [BT [hi [ [k] ] / . FIntin [___ , arg_] ->arg , 

hi [ [maxf] ] / . FIntin [ , arg_] ->arg] , 

maxf =k] (*endlf *) ; 
] (*endFor[k]*) ; 

hi [ [maxf] ] = hi [ [maxf] ] / . FIntin [De___ , _] : > 

FInt [De] ; 

OStep 9(c) *) 

If [Coefficient [eqlist [ [j ] ] ,hl [ [maxf] ] ] ==0 , 
Throw [{eqlist [[j]] ,hl [ [maxf ] ] }] ] Oendif*) ; 
subsrule = {hi [[maxf]] -> 

Collect [ 

(-eqlist [ [j ] ] / . hi [ [maxf] ] ->0) / 
Coefficient [eqlist [ [j] ] ,hl [ [maxf ] ] ] , 

HoldPattern [FInt [ ] ] , Expand [Together] ] } ; 

subslist = Join [subsrule , subslist] ; 

] (*endFor [j]*) ; 

] (*endFor [m2] *) ; 
] (*endFor [ml] *) ; 
] (*endFor[Md2]*) ; 
] (*endFor[Mdl]*) ; 



] (*endFor [1] *) ; 
] (*endFor [Mp] *) ; 
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] OendFor[i2]*) ; 
] OendFor [il] *) ; 

] OendFor [n2]*) ; 
] OendFor [nl] *) ; 
Throw [subsli st] ; 
] OendMod*) ; 
] OendCatch*) ; 
Protect [IBP] ; 
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Appendix B 
Master integrals 



B.l Integrals with up to three external lines 

In this section I give explicit expression for the one-, two- and three-point mas- 
ter integrals, which occur in my calculations. They are calculated in d = 4 — 2e 
dimensions. 

There remains only one nonzero one-point integral: 
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B s2 (x,y) = n 
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Note that ( IB. 21) and (IB. 51) are the leading power of (1B.3I) and ( IB. 61) resp. 
The three-point integrals are: 
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i 1 
(47r) 2 xy 



7T 



ln ?/ 



h In x + In u In y — 

6 y y 2 



— In x ln(x + y — xy) 



+ In y ln(x + y — xy) — 



ln 2 (x + y — xy) 



— In x In # 



+ In ln(x + y - xy) - U 2 x{y X) + Li 2 (y) + Li 2 y( 



x 



-0(0) + 0(e) 



(B.ll) 



B.2 Massive four-point integral 

We consider the following massive four-point integral in d = 4 — 2e dimensions 
(fig. EH : 
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Pi 



k+pi 



P2 




k+pi+ p 2 



Pa/ k+pi+p 2 +P3 



P3 



Figure B.l: Basic one-loop four-point intergral. The massive line, which carries the 
mass m, is indicated by the thick line. 



where 



D 3 



k 2 + it] 

(k + pi) 2 + irj 

{k + pi +P2) 2 +irj 

{k+p! +P2+P3+P4) 2 



m + irj 



(B.13) 



Following [31] we introduce the external masses 



p 2 = m 2 (z = l,2,3,4) 



(B.14) 



and the Mandelstam variables 



S = (Pl+P2) 2 , t=(P2+P3)' 



;b.i5) 



Furthermore we consider only the case, where 



m 2 , = and ml = m 2 . 



(B.16) 



The integral ( IB. 121) can be evaluated using the method of [HI]. This paper gives 
explicit expressions for massless one-loop box integrals. It is however possible to 
extend the single steps of this paper to our case. 



B.3 Massless five-point integral 
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So finally we obtain: 



r(l + e)(4vr/x 2 ) e 



(Aw) 2 m 2 (s — m\) — st + m\m 2 



x 



- (ln(— s — ir]) + ln(m 2 —t — irj) — ln(m 2 — m 2 — irj) — ln(— m\ — ir])) 

+ ln 2 (m 2 — m 2 — irj) + ln 2 (— m 2 — irj) — ln 2 (— s — irj) — ln 2 (m 2 — t — ir]) 

+ ln(m 2 — irj) (ln(— s — irj) + ln(m 2 — t — irj) — ln(m 2 — m 2 — irj) — ln(— m\ — irj)) 

2Li 2 I 1 - m2 -*-.^ - 2Li 2 ( 1 - - mi - * 



— mj — ir] J \ —s — irj 

2Li 2 (1 - (ml - m 2 + zr/)/ m ) + 2Li 2 (l - (mj + i V )f m ) 



2Li 2 (1 - (t - m 2 + ir])f m ) - 2Li 2 (1 - (s + irj)f 



s+t—m\—rn\ 



(B.17) 



where f m = ,, i — m — "t— ?• 

•> mr\m\— s)+st— m^mf 

The case m\ = gives rise to further divergences and has to be considered 
separately: 

.2 n^2^ i r(l + e)(47r/i 2 r 



7 4 (s, t, m 1 = 0, m 3 , m 



x 



3 

2^2 
2tt 2 



(An) 2 s(m 2 - t) 

- ^2 ln(m 2 — t — irj) — — ln(m 2 — irj) + ln(— s — irj) — ln(m 2 — m\ — irj)\ 



+ - ln 2 (m 2 — irj) — ln 2 (m 2 — t — irj) + ln 2 (m 2 — m\ — irj) — ln 2 (— s — irj) 



+ ln(m 2 — irj) (ln(— s — irj) — ln( 



m — m q 



irj)) 



- 2Li 2 1 



m 2 — m\ — ir] 
—s — ir] 



2Li 2 (1 - (m 2 3 - m 2 + zr])f m ) 



2Li 2 (1 - (t - m 2 + ir])f m ) - 2Li 2 (1 - (s + mj)/' 



(B.18) 



where f r 



s(t—m 2 ) ' 



B.3 Massless five-point integral 

We consider the following massless five-point integral: 

d d k 1 



E = ]X 



2c 



(2vr) d k 2 (k + xp) 2 (k + p) 2 (k + l) 2 (k + l- yq) 2 ' 



(B.19) 
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B. Master integrals 



Despite the fact that there occur only three linearly independent momenta in E , 
(IB. 191) cannot be decomposed into partial fractions. However using a standard 
Feynman parametrisation we end up with integrals which can be calculated by 
computer algebra systems. Now the exact result is rather involved so I give just 
the leading power. Higher powers can be obtained by using the methods of section 



% r(l + e)(4vr/i 2 ) 



2 -4 + 21nf + 21n# + 2i7r 

-- + " 

e l e 

- 21n£ln# + 41n£ + 41n# + 7r 2 



In 2 f - In 2 9 



o — 



(4tt) 2 xydi 2 



+ 2in (2-21nf -21n0) 



(B.20) 



Appendix C 
Matching of 



In this thesis we consider a pure QCD calculation. This calculation is given in terms 
of Jb/^b- In order to compare our calculation to HQET or SCET results we have 
to match the expression /b/Ab onto HQET. We use the definition of the .B-meson 
wave function of [35], which is defined in the case of HQET by the heavy quark field 
of h v . So we define 

1 r°° 

ifBmg^iu) = — J dte^(0\q(z)[. . .]ji +l5 b(0)\B) z - <z±=0 (C.l) 
in the case of QCD and 

i r°° 

ifBm B <f>* Q * T (u) = — J dte^(0\q(z)[...]ji +l5 h v (0)\B) z -, z± =o (C.2) 

in the case of HQET. Here the integration goes over t = v ■ z where v is the four- 
velocity of the 5-meson. We define the 5-meson decay constant by 

if B m B = (0\q(0)[--]^5K0)\B) (C.3) 

and analogously the HQET decay constant, which depends on the renormalisation 
scale (j,, by 

if B (»)m B = (0|g(0)[. . ]i l5 h v (0)\B). (C.4) 
The matching coefficient is defined by 

Jo u J UJ 

such that we get writing the /x-dependence explicitly 

fB -CxMtM 1 -- (C6) 



A QCD (//) — A HQET (/i) - 

We get C\ B to 0(a s ) by calculating the convolution integrals occurring in ( 1C.5I) 
in both QCD and HQET up to 0{a s ). The corresponding diagrams are shown in 
fig. IC.ll As in section I3T31 we can use the wave functions flClH . ( 1C.2I) defined by free 
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Figure C.l: NLO contributions to X B . The double line stands for the 6-quark field. 

quark states instead of hadronic states. We assign to the 6-quark the momentum 
v(rrib — Co) (—vCo resp.) in the case of pure QCD (HQET resp.) and vCo to the soft 
constituent quark, where v is the four velocity of the S-meson. At tree level we get 
for both QCD and HQET the same wave function: 

if B <f> { +\co) = N c 5(uj - Co) qfi +l5 ^ (C.7) 

where the spinor \l/ fulfils the condition pftp = and our convolution integral is 

if B r du = _l ^ 75 *. (c.8) 

Jo oj N c u 

At NLO only the first diagram in fig. IC.ll needs to be considered as the other two 
are in leading power identical for QCD and HQET. The following expressions are 
given in the MS scheme, i.e. we redefine fi 2 — > f^ 2 ^-- For the diagram in fig. IC.lf a) 
we get in QCD 

^C F N c lqi^ f 2 + 21n ^ +41nA + 4 -^-21n 2 ^ + 41n^lnA 

47r uj \ e rrib 6 nib mj m& 

(C.9) 

and in HQET 

a 1 ( 1 21n^ n u tt 2 \ 

Air uj \ € e fi 4 I 

The wave function renormalisation constants of the heavy quark field are given in 
the onshell scheme for the QCD 6-field: 

Z\ = l + ^CV (-—- — -3 In -^-2^) (C.ll) 
26 4tt F \ 2e e m m b J 1 ; 

and for the HQET field h v : 
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The renormalisation of the g-field drops out in the matching. Diagrammatically the 
matching equation (10. 5|) reads: 



72 

L 2b 



/ 








\ 


V 




+ 




/ 



QCD 



/ 



C\ B Z 2h v 



HQET 



(C.13) 



Finally we obtain 



ft. 



47T 



C Xg (fj) = 1 + ^C F ( 2 In 2 A + In + 2 + — 



7T 



12 



(C.14) 



where we have renormalised the UV-divergences in the MS-scheme. 



104 C. Matching of \ B 



Bibliography 



[1] L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 

[2] A. J. Buras, M. E. Lautenbacher and G. Ostermaier, Phys. Rev. D50, 3433 
(1994), [hep-ph/9403384]. 

[3] CKMfitter Group, J. Charles et al, Eur. Phys. J. C41, 1 (2005), [hep- 
ph/0406184], updates at http://ckmfitter.in2p3.fr. 

[4] BELLE, K. Abe et al, Phys. Rev. D71, 072003 (2005), [hep-ex/0408111]. 

[5] BABAR, B. Aubert et al, Phys. Rev. D73, 012004 (2006), [hep-ex/0507054]. 

[6] BABAR, e. Harrison, P. F. and e. Quinn, Helen R., Papers from Workshop on 
Physics at an Asymmetric B Factory (BaBar Collaboration Meeting), Rome, 
Italy, 11-14 Nov 1996, Princeton, NJ, 17-20 Mar 1997, Orsay, France, 16-19 Jun 
1997 and Pasadena, CA, 22-24 Sep 1997. 

[7] Heavy Flavor Averaging Group (HFAG), E. Barberio et al, [hep-ex/0603003], 
updates at http://www.slac.stanford.edu/xorg/hfag. 

[8] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Phys. Rev. Lett. 
83, 1914 (1999), [hep-ph/9905312]. 

[9] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Nucl. Phys. B591, 
313 (2000), [hep-ph/0006124]. 

[10] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Nucl. Phys. B606, 
245 (2001), [hep-ph/0104110]. 

[11] G. Bell, [arXiv:0705.3133 [hep-ph]]. 

[12] G. Bell, [arXiv:0705.3127 [hep-ph]]. 

[13] M. Beneke and S. Jager, Nucl. Phys. B751, 160 (2006), [hep-ph/0512351]. 
[14] N. Kivel, [hep-ph/0608291]. 

[15] C. W. Bauer, S. Fleming, D. Pirjol and I. W. Stewart, Phys. Rev. D63, 114020 
(2001), [hep-ph/0011336]. 



106 



BIBLIOGRAPHY 



[16] C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. D65, 054022 (2002), 
[hep-ph/0109045]. 

[17] M. Beneke, A. P. Chapovsky, M. Diehl and T. Feldmann, Nucl. Phys. B643, 
431 (2002), [hep-ph/0206152]. 

[18] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 
(1996), [hep-ph/95 12380]. 

[19] M. Gronau, O. F. Hernandez, D. London and J. L. Rosner, Phys. Rev. D50, 
4529 (1994), [hep-ph/9404283]. 

[20] K. G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B192, 159 (1981). 

[21] F. V. Tkachov, Phys. Lett. B100, 65 (1981). 

[22] S. Laporta, Int. J. Mod. Phys. A15, 5087 (2000), [hep-ph/0102033]. 

[23] G. Passarino and M. J. G. Veltman, Nucl. Phys. B160, 151 (1979). 

[24] T. Gehrmann and E. Remiddi, Nucl. Phys. B580, 485 (2000), [hep- 
ph/9912329]. 

[25] E. Remiddi, Nuovo Cim. A110, 1435 (1997), [hep-th/9711188]. 

[26] M. Caffo, H. Czyz, S. Laporta and E. Remiddi, Nuovo Cim. Alll, 365 (1998), 
[hep-th/9805118]. 

[27] S. G. Gorishnii, Nucl. Phys. B319, 633 (1989). 

[28] M. Beneke and V. A. Smirnov, Nucl. Phys. B522, 321 (1998), [hep- 
ph/9711391]. 

[29] V. A. Smirnov, Commun. Math. Phys. 134, 109 (1990). 

[30] V. A. Smirnov, Springer Tracts Mod. Phys. 177, 1 (2002). 

[31] G. Duplancic and B. Nizic, Eur. Phys. J. C20, 357 (2001), [hep-ph/0006249]. 

[32] S. Herrlich and U. Nierste, Nucl. Phys. B455, 39 (1995), [hep-ph/9412375]. 

[33] A. F. Falk, B. Grinstein and M. E. Luke, Nucl. Phys. B357, 185 (1991). 

[34] G. P. Lepage and S. J. Brodsky, Phys. Rev. D22, 2157 (1980). 

[35] B. O. Lange and M. Neubert, Phys. Rev. Lett. 91, 102001 (2003), [hep- 
ph/0303082]. 

[36] M. Jamin and B. O. Lange, Phys. Rev. D65, 056005 (2002), [hep-ph/0108135]. 
[37] A. Khodjamirian, T. Mannel and N. Offen, [hep-ph/0611193]. 



BIBLIOGRAPHY 



107 



[38] A. Abada et ai, Nucl. Phys. B619, 565 (2001), [hep-lat/0011065]. 

[39] E. Dalgic et ai, Phys. Rev. D73, 074502 (2006), [hep-lat/0601021]. 

[40] A. Khodjamirian, R. Riickl, S. Weinzierl, C. W. Winhart and O. I. Yakovlev, 
Phys. Rev. D62, 114002 (2000), [hep-ph/0001297]. 

[41] M. G6ckeler et ai, Nucl. Phys. Proc. Suppl. 161, 69 (2006), [hep-lat/0510089]. 

[42] V. M. Braun et ai, Phys. Rev. D74, 074501 (2006), [hep-lat/0606012]. 

[43] Particle Data Group, W. M. Yao et al, J. Phys. G33, 1 (2006), updates at 
http : //pdg . lbl . gov. 

[44] A. Khodjamirian, T. Mannel and N. Offen, Phys. Lett. B620, 52 (2005), [hep- 
ph/0504091]. 

[45] V. M. Braun, D. Y. Ivanov and G. P. Korchemsky, Phys. Rev. D69, 034014 
(2004), [hep-ph/0309330]. 

[46] A. G. Grozin and M. Neubert, Phys. Rev. D55, 272 (1997), [hep-ph/9607366]. 

[47] S. J. Lee and M. Neubert, Phys. Rev. D72, 094028 (2005), [hep-ph/0509350]. 

[48] M. Beneke and M. Neubert, Nucl. Phys. B675, 333 (2003), [hep-ph/0308039]. 

[49] M. Beneke and S. Jager, Nucl. Phys. B768, 51 (2007), [hep-ph/0610322]. 

[50] M. Jamin and M. E. Lautenbacher, Comput. Phys. Commun. 74, 265 (1993). 



108 



Acknowledgements 



There are many people who contributed to the present thesis. First of all I would 
like to thank my advisor Prof. Gerhard Buchalla. He made my PhD studies possible 
by supplying me with a scholarship of the "Graduiertenkolleg Particle Physics at the 
Energy Frontier of New Phenomena" . I also want to thank him for many inspiring 
discussions. 

I am particularly indebted to Guido Bell for giving me useful hints and for stim- 
ulating discussions, where I learned much about technical details in the calculations 
of Feynman integrals. I would also like to thank my room mate Matthaus Bartsch 
and Sebastian Jager for many instructive and helpful discussions. 

Furthermore I want to thank Matthaus Bartsch, Guido Bell, Gerhard Buchalla 
and Sebastian Jager for proofreading the drafts and comments on the manuscript. 

Finally I want to thank my parents, whose aid made my physics studies possible. 



